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Abstract 

This  thesis  examines  the  ground  and  excited  electronic  states  of  the  uranyl 
(UO2  )  and  uranate  (UO4  ')  ions  using  Hartree-Fock  self-consistent  field  (HF  SCF), 
multi-configuration  self-consistent  field  (MCSCF),  and  multi-reference  single  and  double 
excitation  configuration  interaction  (MR-CISD)  methods.  The  MR-CISD  calculation 
included  spin-orbit  operators.  Molecular  geometries  were  obtained  from  self-consistent 
field  (SCF),  second-order  perturbation  theory  (MP2),  and  density  functional  theory 
(DFT)  geometry  optimizations  using  the  NWChem  4.01  massively  parallel  ab  initio 
software  package.  COLUMBUS  version  5.8. 1  was  used  to  perform  in-depth  analysis  on 
the  HF  SCF,  MCSCF,  and  MR-CISD  potential  energy  surfaces. 

Excited  state  calculations  for  the  uranyl  ion  were  performed  using  both  a  large- 
and  small  core  relativistic  effective  core  potential  (RECP)  in  order  to  calibrate  the 
method.  This  calibration  included  comparison  to  previous  theoretical  and  experimental 
work  on  the  uranyl  ion.  Uranate  excited  states  were  performed  using  the  small-core 
RECP  as  well  as  the  methodology  developed  using  the  uranyl  ion. 
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THEORETICAL  COMPARISON  OF  THE  EXCITED  ELECTRONIC  STATES  OF  THE 


LINEAR  URANYL  (U022+)  AND  TETRAHEDRAL  URANATE  (U042')  IONS  USING 
RELATIVISTIC  COMPUTATIONAL  METHODS 


I.  Introduction 


The  chemical  properties  of  uranium  and  plutonium  oxides  are  critically  important 
to  nuclear  applications.  Of  particular  importance  in  the  chemistry  of  these  actinide 
compounds  is  the  oxidation  state.  Uranium,  like  all  the  actinides,  can  possess  a  wide 
range  of  oxidation  states,  ranging  from  +3  to  +6,  due  to  chemical  activation  of  the 
uranium  5f  orbitals  via  relativistic  effects.  As  a  result,  the  oxidation  state  of  uranium  can 
be  influenced  by  its  local  chemical  environment,  which  in  turn  influences  the  geometry 
of  the  uranium  oxide  compounds.  Uranium  oxidation  state  plays  a  very  important  role 
throughout  the  nuclear  fuel  cycle,  and  it  plays  a  critical  role  in  the  mobility  of  uranium  in 
the  environment.  This  oxidation  state  can  be  inferred  through  spectroscopic 
measurements,  providing  a  simple  and  inexpensive  tool  for  use  in  such  areas  as  nuclear 
forensics  and  environmental  monitoring. 

Additionally,  the  stockpile  stewardship1  program  demands  a  thorough 

understanding  of  the  processes  by  which  uranium  and  plutonium  components  age,  as  well 

as  the  effect  this  aging  has  on  the  reliability  and  perfonnance  of  nuclear  weapons.  A 

cornerstone  of  the  stockpile  stewardship  program  is  theoretical  modeling  and  simulation 

1  Stockpile  stewardship  refers  to  the  substantial  effort  undertaken  by  the  U.S.  Department  of  Energy  to 
maintain  and  certify  the  U.S.  nuclear  weapon  arsenal  without  resorting  to  underground  nuclear  testing— 
http://www.doe.gov. 
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of  the  basic  physics  and  chemistry  involved  in  the  design,  manufacture,  maintenance,  and 
operation  of  a  nuclear  weapon.  Again,  non-invasive  electronic  spectroscopic  methods 
can  be  used  to  diagnose  the  extent  of  nuclear  weapon  component  aging,  based  on  the 
relationship  between  uranium  oxidation  state  and  its  local  chemical  environment. 

Uranium  oxidation  is  of  particular  interest.  Oxygen  and  uranium  readily  react, 
forming  a  wide  range  of  complex  oxides  over  a  range  of  temperatures  and  pressures 
(Wanner,  1992).  The  uranyl  ion,  U02  ,  is  an  unusually  stable  oxide  of  uranium,  and  it  is 
present  in  a  majority  of  complex  uranium(VI)  oxides  (Pyykko,  1998:  3787-3788;  Zhang 
1999:6880). 

While  there  is  a  large  amount  of  experimental  data  on  the  various  properties  of 
uranium  and  plutonium  (Katz  et  al,  1986;  Wanner,  1992),  theoretical  understanding  of 
the  spectra  of  these  elements  has  progressed  slowly.  Ab  initio 2  quantum  mechanical 
theoretical  techniques  have  made  great  strides  in  understanding  of  molecules  consisting 
of  lighter  elements,  and  computational  methods  have  been  quite  successful  in  predicting 
thennodynamic  and  spectroscopic  properties  of  these  compounds.  Unfortunately,  such 
progress  in  the  actinide  compounds  has  progressed  more  slowly,  for  two  main  reasons. 

The  first  difficulty  is  the  sheer  number  of  electrons  to  deal  with  in  actinide 
compounds.  Common  uranium  oxide  compounds  such  as  U02  have  108  electrons,  while 
more  complex  oxides  such  as  U3O8  have  over  300  electrons.  Accurately  treating  such 
large  numbers  of  electrons  becomes  computationally  intensive,  and  it  has  only  been  in 
the  last  decade  that  such  molecules  can  be  treated  with  the  accuracy  necessary  to 
compare  theoretical  and  experimental  electronic  spectra.  A  second  difficulty  is  the  fact 
that  relativistic  effects  must  be  accounted,  not  as  perturbations  to,  but  on  an  equal  footing 
with  electron  correlation  in  these  heavy  molecular  systems  for  even  moderate  accuracy. 
This  is  in  stark  contrast  to  lighter  molecules  where  relativistic  effects  can  be  neglected  in 

2  Latin  for  "from  the  beginning",  Merriam-Webster's  Collegiate  Dictionary,  2002,  http://www.webster.com 
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all  but  high-precision  theoretical  calculations  (Pyykko,  1998:  563-594;  Balasubramanian, 
1997:  1-27). 

With  the  development  of  faster  computers,  especially  massively  parallel  computer 
systems,  as  well  as  quantum  chemistry  software  codes  designed  to  take  advantage  of 
these  computer  architectures,  there  has  been  an  increasing  number  of  theoretical 
investigations  of  actinide  compounds.  However,  the  theoretical  calculation  of  excited 
electronic  states  for  actinide  compounds  is  a  difficult  problem,  and  there  are  relatively 
few  in-depth  studies  of  the  excited  states  of  uranium  oxides  necessary  for  understanding 
the  electronic  spectra. 

This  research  focuses  on  two  uranium  oxide  ions  in  particular:  uranyl  (UO2  ) 
and  uranate  (U04“‘)  ions.  Starting  from  optimized,  theoretical,  gas-phase  molecular 
geometries,  electronic  spectra  calculations  from  single  and  double  excitations  with  spin- 
orbit  coupling  included  were  computed  and  compared  with  experimental  and  other 
theoretical  results.  The  calculations  involving  uranyl  were  used  to  calibrate  and  validate 
the  method,  while  those  involving  U04  'were  an  attempt  to  begin  understanding  the 
influence  of  the  local  oxygen  coordination  on  the  electronic  spectra  of  uranium  oxides. 
Particular  interest  was  paid  to  the  first  excited  states  of  both. 

The  theory  relevant  to  calculations  of  the  electronic  spectra  of  uranium  oxides  is 
laid  out  in  chapter  two.  Next,  the  hardware  and  software  resources  used  in  this  research, 
as  well  as  the  methodology  behind  the  study  of  uranyl  and  uranate  electronic  spectra  is 
described  in  chapter  three.  Results  and  discussion  of  the  results  is  included  in  chapter 
four,  followed  by  conclusions  drawn  from  this  research  and  recommendations  for  further 
research  in  chapter  five. 
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II.  Theory 


Ab  initio  quantum  mechanical  theoretical  techniques  have  been  applied  to 
molecules  containing  light  elements  with  increasing  success  in  the  past  several  decades. 
Application  of  these  theoretical  techniques  to  light-atom  molecules,  especially  organic 
molecules,  has  yielded  vast  insight  into  the  properties  of  these  molecules  such  as 
molecular  geometries  for  ground  and  transition  states,  electron  affinities,  ionization 
potentials  and  more.  Advances  in  computing  power,  coupled  with  quantum  chemistry 
software  designed  specifically  to  take  maximum  advantage  of  these  computers  has 
increased  the  applicability  to  larger  molecules.  Such  calculations  have  become  an 
indispensable  tool  to  theoreticians  and  experimentalists  alike. 

Complications  arise  when  applying  theory  to  molecules  containing  heavy 
elements,  especially  actinide  molecules.  The  two  most  difficult  complications  to  the 
theoretical  treatment  of  heavy-element  molecules  are  increased  electron  correlation  and 
relativistic  effects. 

First,  these  heavy-element  molecules  contain  a  large  number  of  electrons  whose 
motions  are  coupled  through  electrostatic  and  quantum  mechanical  interactions.  Electron 
correlation  effects  can  contribute  roughly  1  eV  (23  kcal/mol)  to  the  total  electronic 
energy  per  electron  pair  (Raghavachari,  1996:  12960).  Using  this  rule  of  thumb, 
electron  correlation  accounts  for  46  eV  of  the  total  electronic  energy  in  the  uranium 
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atom.  An  accurate  treatment  of  electronic  correlation  is  critical  in  order  to  perform 
meaningful  comparisons  between  theoretical  and  experimental  spectra. 

A  second  complicating  factor  in  the  theoretical  treatment  of  heavy-element 
molecules  is  the  increasing  importance  relativistic  effects  play  in  the  accurate  description 
of  ground  and  excited  electronic  states  with  increasing  atomic  number,  Z.  Several 
reviews  examine  relevant  chemical  effects  due  to  relativistic  quantum  mechanical 
treatments  (Pyykko,  1988:  563-594;  Pepper  et  al,  1991:  719-741;  Kaltsoyannis,  1997:  1- 
11). 


Relativistic  Effects  in  Chemistry 

There  are  three  main  relativistic  effects  in  atomic  and  molecular  chemistry,  all  of 
which  are  roughly  the  same  magnitude,  and  they  approximately  scale  as  Z4  (Pyykko, 

1988,  564).  The  first  main  relativistic  effect  is  considered  a  direct  relativistic  effect,  and 
it  consists  of  a  radial  contraction  of  atomic  orbitals,  along  with  a  lowering  of  the  energy 
level  of  the  electronic  state.  This  effect  is  due  primarily  to  the  relativistic  mass  increase 
as  electron  velocities  become  appreciable  fractions  of  the  speed  of  light.  Simple 
replacement  of  the  relativistic  mass  expression  for  the  electron  in  the  Bohr  radius  formula 
yields 


4  ns0h2^Jl-(^)2 


(1) 


Here,  h  is  Planck’s  constant  divided  by  2jt,  s0  is  the  pennittivity  of  free  space,  e 
is  the  electron  charge,  and  m0  is  the  electron  mass.  As  electron  speeds,  v,  approach  the 
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speed  of  light,  c,  the  Bohr  radius,  ao,  shrinks.  Electron  orbitals  with  high  densities  near 
the  nucleus  experience  the  largest  contractions,  where  electron  speeds  are  largest.  For 
electrons  in  the  hydrogenic  Is  shell,  the  average  fraction  of  the  speed  of  light  is  given  by 
Z 

(Pyykko,  1988:  563).  For  uranium,  this  is  0.67c,  yielding  a  Is  orbital  radial 

contraction  of  roughly  26%.  Ah  atomic  orbitals  have  some  density  near  the  nucleus, 
therefore,  ah  atomic  orbitals  experience  some  contraction.  However,  the  inner  s-  and  p- 
orbitals  nearest  the  nucleus  experience  the  most  contraction  (Pyykko,  1988:  563).  In 
light-element  molecules,  this  orbital  contraction  is  small  and  negligible  in  ah  but  the 
highest  precision  calculations,  but  the  effect  becomes  dramatic  in  actinide  elements  such 
as  uranium. 

The  second  relativistic  effect  is  considered  to  be  an  indirect  effect,  and  it  consists 
of  a  radial  expansion  and  increase  in  the  electronic  energy  levels  of  outer  atomic  orbitals. 
This  is  due  to  more  effective  nuclear  charge  screening  by  the  inner,  contracted  electrons, 
reducing  the  effective  nuclear  charge  experienced  by  the  outer  electrons.  Additionally, 
relativistic  contraction  of  the  inner  s-  and  p-  electron  shells  increase  the  electron  density 
near  the  nucleus,  crowding  out  the  outer  d-  and  f-  electron  shells.  This  is  due  to  the  fact 
that  there  is  a  decrease  in  electron  density  near  the  nucleus  for  orbitals  with  increasing 
orbital  angular  momentum.  Thus,  the  direct  orbital  contraction  competes  with  the 
indirect  orbital  expansion.  In  general,  the  result  of  this  interplay  between  relativistic 
effects  is  to  contract  and  stabilize  s-  and  p-  atomic  orbitals,  while  d-  and  f-  orbitals 
expand  and  destabilize  in  energy.  The  orbital  expansion  and  contraction  can  affect  bond 
lengths  (Pyykko,  1988:  571)  and  force  constants  (Pyykko,  1988:  580),  which  in  turn 
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affect  molecular  vibrational  frequencies.  These  relativistic  effects  obviously  affect  the 
observed  spectra  of  heavy-element  molecules,  but  not  as  much  as  the  splitting  of  states 
caused  by  the  third  major  relativistic  effect,  spin-orbit  coupling. 

Intrinsic  electron  spin  is  a  natural  result  of  a  Lorentz-covariant  description  of  the 
quantum  mechanically  wave  equation  (Balasubramanian,  1997:  76-78).  This  spin 
angular  momentum  couples  with  the  electron  orbital  angular  momentum,  lifting 
degeneracy  in  atomic  orbitals  with  angular  momentum.  Thus,  the  three  degenerate  p- 

orbitals  in  non-relativistic  theory  split  into  one  p,  orbital  and  two  degenerate  p  . 

2  2 

orbitals.  Of  the  three  effects,  spin-orbit  coupling  has  the  largest  impact  in  atomic  and 
molecular  spectra,  even  for  low-Z  atoms  and  molecules.  For  light  atoms,  a  perturbative 
treatment  of  spin-orbit  coupling  known  as  Russell-Sanders  coupling  or  L-S  coupling 
often  yields  sufficient  accuracy  for  electronic  transition  energies.  This  coupling  scheme 
treats  magnetic  spin-orbit  coupling  as  a  small  perturbation  to  the  electron-electron 
electrostatic  interaction.  Orbital  angular  momentum  and  spin  angular  momentum  are  still 
“nearly”  good  quantum  numbers  in  this  coupling  scheme,  and  both  L  and  S  commute 
with  the  Hamiltonian  in  Russell-Sanders  coupling  scheme.  Atomic  states  are  described 
by  term  symbols 

25+1  j 
^ J 

S  =  0,1,2,... 

L  =  S,P,D,F ,... 

J  =  L  +  S 

with  S  equal  to  the  total  spin  multiplicity,  and  L  is  the  total  orbital  angular  momentum  (0, 
1,  2,  ..).  Traditional  spectroscopic  notation  is  used  for  the  total  orbital  angular 
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momentum,  with  S  representing  zero  total  orbital  angular  momentum,  P  representing  one 
unit  of  orbital  angular  momentum  and  so  on.  J  is  the  total  angular  momentum  of  the 
electron,  given  by  the  sum  of  orbital  and  spin  angular  momenta.  Examples  of  Russell- 
Sanders  tenn  symbols  include  3P2,  1  So,  and  3D4  (Gerloch,  1986:  69-74). 

On  the  other  end  of  the  perturbation  spectrum,  more  appropriate  for  very  heavy 
atoms,  the  electron-electron  electrostatic  interaction  is  treated  as  a  perturbation  to  the 
magnetic  spin-orbit  coupling.  This  coupling  scheme  is  known  as  j-j  coupling.  In  this 
coupling  scheme,  neither  L  nor  S  commute  with  the  Hamiltonian.  However,  the  total 
angular  momentum,  J,  still  commutes  with  the  atomic  Hamiltonian,  and  hence,  is  a  good 
quantum  number.  The  term  symbol  for  j-j  coupling  is  given  by  the  J  value  for  the  state 
(Gerloch,  1986:  74-76). 

Figure  1  contains  a  schematic  representation  of  the  two  spin-orbit  coupling 
extremes  (Gerloch,  1986:  61).  The  horizontal  springs  represent  electrostatic  coupling 
between  the  electrons,  while  the  vertical  springs  represent  magnetic  coupling  between  the 
electron  intrinsic  magnetic  moments.  In  Russell-Sanders  coupling,  the  electron  orbital 
angular  momenta  couple  strongly,  as  do  each  electron’s  spin  angular  momenta.  These 
total  orbital  and  spin  angular  momenta  then  couple  weakly.  The  opposite  is  true  of  j-j 
coupling.  In  j-j  coupling,  each  electron’s  orbital  and  spin  angular  momenta  couple 
strongly,  and  this  individual  total  angular  momentum  couples  weakly  with  the  other 


electrons  total  angular  momenta. 


Coulomb 


coupling 

scheme 

Figure  1.  Schematic  Representation  of  the  Assumptions  in  Russell-Sanders  and  j-j  Spin- 
Orbit  Coupling  Schemes  (Gerloch,  1986:  61) 

Most  elements  on  the  periodic  table  fall  between  these  two  perturbation  extremes, 
and  so  intennediate  coupling  is  more  appropriate  than  either  perturbative  treatment. 
Intermediate  coupling  is  not  a  separate  coupling  scheme,  but  occurs  as  deviations  from 
the  separate  perturbative  treatments  given  by  L-S  and  j-j  coupling  (Gerloch,  1986:  77). 
Figure  2  illustrates  the  effect  of  Russell-Sanders,  intermediate,  and  j-j  spin-orbit  coupling 
on  a  d2  electronic  configuration.  Figure  2  shows  the  effect  of  spin-orbit  coupling  on  an 
atomic  electronic  state  with  two  electrons  in  the  d-shell.  The  left-hand  side  shows  the 
tenn  symbols  that  arise  due  to  Russell-Sanders  coupling,  while  the  right-hand  side  shows 
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the  effect  of  j-j  coupling  on  the  same  electronic  configuration.  Here,  the  importance  of 
spin-orbit  coupling  to  electronic  spectroscopy  is  evident.  Without  taking  into  account 
spin-orbit  coupling,  both  the  number  of  states  and  their  relative  ordering  will  be 
incorrect. 


j 


% 

R-S  coupling  Intermediate  coupling  J  -  j  coupling 
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Figure  2.  Correlation  Diagram  of  the  Various  States  Arising  From  a  d“  Electronic 
Configuration  Using  Both  Russell-Sanders  and  j-j  Spin-Orbit  Coupling  (Gerloch,  1986: 

78) 

Relativity  also  affects  the  symmetry  of  molecules,  because  of  electron  spin. 
Under  the  assumption  that  the  total  electronic  wave  function  can  separated  (the  product 
of  a  spatial  and  spin  wave  functions),  each  wave  function  may  possess  separate 
symmetry,  and  the  total,  observable  state  symmetry  is  given  by  the  product  of  the  spatial 
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and  spin  symmetries.  For  example,  singlet  spin  states  are  completely  symmetric,  while 
triplet  spin  states  transfonn  like  the  components  of  the  angular  momentum  operator. 

Thus,  a  completely  symmetric  spatial  wave  function  multiplied  by  a  triplet  spin  wave 
function  will  not  be  totally  symmetric.  For  systems  with  a  spin-orbit  Hamiltonian,  the 
symmetry  point  groups  can  have  twice  the  number  of  symmetry  operations,  and  are 
called  double  point  groups.  This  doubling  of  the  order  of  the  symmetry  point  groups  is 
due  to  the  introduction  of  \  -integral  angular  momentum  values.  Systems  possessing  an 

even  number  of  electrons  obey  Bose-Einstein  (bosons)  statistics,  and  the  total  wave 
function  of  bosonic  systems  is  symmetric  with  respect  to  rotations  by  2n.  Systems 
possessing  an  odd  number  of  electrons  obey  Fermi-Dirac  statistics  (fermions),  and 
fermionic  wave  functions  change  sign  up  the  exchange  of  two  particles.  This  exchange  is 
equivalent  to  a  rotation  by  2jt,  and  so  a  rotation  of  4ji  returns  a  fermionic  system  to  its 
original  state.  While  bosonic  systems  transform  according  to  the  irreducible 
representations  of  the  single  point  groups,  the  rotation  by  2%  is  a  new  symmetry  operation 
for  fermionic  systems,  doubling  the  order  of  the  symmetry  point  group.  For  example, 
rotations  of  a  closed-shell  molecule,  such  as  uranyl  (UO2  )  transforms  according  to  the 
nonnal  irreducible  representations  of  the  Dooh  point  group.  Rotating  the  molecule  by  2n 
leaves  the  molecule  (wave  function)  unchanged.  However,  for  an  open-shell  molecule, 
such  as  UO2,  such  is  not  the  case.  Such  molecules  transform  according  to  the  extra 
irreducible  representations  generated  by  a  rotation  of  2n.  Rotating  the  UO2  molecule  by 
27i  introduces  a  phase  factor  into  the  total  electronic  wave  function.  A  rotation  by  4tt  in 
this  case  returns  the  molecule  (wave  function)  to  its  original  configuration. 
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These  three  main  effects,  orbital  contraction  and  energy  stabilization,  orbital 
expansion  and  energy  destabilization,  and  spin-orbit  coupling,  along  with  the 
consequential  double  group  symmetry  constitute  the  chemically  relevant  relativistic 
effects  in  atoms  and  molecules.  The  most  important,  from  a  spectroscopic  standpoint  is 
spin-orbit  coupling,  even  in  the  spectra  of  the  lightest  elements.  A  quantum  mechanical 
treatment  of  the  electron  must  account  for  the  intrinsic  magnetic  moment  of  the  electron, 
and  the  Dirac  equation  accomplishes  this  quite  elegantly. 


The  Dirac  Equation 

Relativity  has  played  a  role  in  quantum  mechanical  systems  since  the  inception  of 
the  theory.  Attempts  at  finding  a  Lorentz  invariant  form  for  Schrodinger’s  equation  led 
to  two  Lorentz-covariant  equations:  the  Klein-Gordon  equation,  and  the  Dirac  equation. 
Schrodinger’s  equation,  a  non-relativistic  quantum  mechanical  wave  equation,  is  given 
by 

IW  =  EXY  ,  (2) 

where  the  Hamiltonian,  H,  and  energy,  E,  operators  are  given  by 

H{f,t)  =  E^A  +  V{7,t),  (3) 

2m 


and 


m= 


h  d 
i  St 


The  momentum  operator,  p,  is  defined  by 

p(r,t)  =  -V. 

i 


(4) 


(5) 
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Here,  m  is  the  electron  mass,  i  is  the  imaginary  number,  fi  is  Planck’s  constant  divided 
by  2k  ,  V(r,t)  is  the  potential  energy  operator,  and  p(r,t )  is  the  electron  momentum. 

'P  is  the  electronic  wave  function. 

This  equation,  because  of  the  non-equivalent  treatment  of  the  spatial  and  temporal 
variables,  is  not  Lorentz  invariant,  and  therefore  is  limited  to  non-relativistic  phenomena. 
Early  attempts  at  making  a  Lorentz-covariant  equation  began  by  quantizing  the  Lorentz- 
covariant  relativistic  energy  expression 

E1  =  p2c2  +m20c\  (6) 

Again,  p  is  the  electron  momentum,  c  is  the  speed  of  light,  and  m0  is  the  electron 
rest-mass,  and  E  is  the  electron  energy. 

Replacing  the  energy  and  momentum  expressions  with  their  quantized 
counterparts  leads  to  the  Klein-Gordon  wave  equation  for  a  free  particle 
(Balasubramanian,  1997:  99-101;  Messiah,  1999:  884-888;  Bjorken  and  Drell,  1964:  4- 
6,198-206) 

-n2^-  =  -n2c2V2x¥ +  myx¥ .  (7) 

dt2 

While  this  scalar  wave  function  is  Lorentz-covariant,  it  has  several  undesirable 
properties,  making  it  unacceptable  as  a  wave  function  for  the  electron.  First,  the 
probability  density  associated  with  it  is  not  positive  definite,  resulting  in  possible 
negative  probability  densities.  Additionally,  both  positive  and  negative  energy  solutions 
to  this  equation  exist,  complicating  early  interpretation  of  the  wave  function.  The  fact 
that  the  probability  density  is  not  positive  definite  makes  this  equation  a  poor  choice  for 
an  electronic  wave  function;  however,  the  Klein-Gordon  turns  out  to  be  a  valid 
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relativistic  wave  equation  for  spin- free  fields,  such  as  pi  mesons  (Messiah,  1999:  888; 
Balasubramanian,  1997:  108). 

Dirac  took  a  different  approach  in  formulating  a  Lorentz-covariant  equation  for  a 
free  electron  (Dirac,  1928:  610-624;  Dirac,  1928:  351-361;  Balasubramanian,  1997: 
110-119).  He  began  with  the  same  Lorentz-covariant  expression  for  the  energy  of  a  free 
particle  as  the  Klein-Gordon  equation, 

p 2  22.  24  /o\ 

h  =  p  c  +m0c  .  (8) 

Taking  the  square  root  yields  the  Dirac  Hamiltonian, 

.  (9) 

Quantizing  this  expression  by  the  usual  substitutions  for  the  energy  and 
momentum  operators  yields  a  Hamilton  that  involves  a  first-order  time  derivative. 
However,  the  square  root  in  the  operator  makes  application  problematic  and  hopelessly 
complicated.  Dirac  circumvented  this  problem  by  introducing  a  new  degree  of  freedom 
into  the  relativistic  Hamiltonian,  effectively  completing  the  square.  This  yielded  a  more 
tractable  Hamiltonian  operator 

H Dirac  =  C(«1  •  Pi  +  a2  ■  Pi  +  a3  •  Pi)  +  •  (10) 

The  Dirac  equation, 

ifi  dKV 

- —  =  c(al  ■pl+a2-p2  +  a3  ■  p3)W  +  pm0c2xV ,  (11) 

c  ot 

results  from  this  Hamiltonian  (Kellogg,  1997:  4-6). 

Requiring  solutions  to  this  equation  to  simultaneously  satisfy  the  Klein-Gordon 
equation  places  restrictions  on  the  components  of  the  a,  and  ft  matrices: 
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aiaJ  +  a  jai  =  28 tj , 


(12) 


P2  =  1,  (13) 

and 

akP  +  pak=  0.  (14) 

In  order  to  satisfy  these  restrictions  both  ai  and  P  must  be  at  least  four-by-four 
matrices,  which  operate  on  a  four-component  vector  wave  function.  The  ai  and  P 
matrices  are  defined  as 

0  0  0  1 
0  0  10 
0  10  0 
10  0  0 

0  0  0  -/ 

0  0/0 
0-/0  0 
/  0  0  0 

0  0  10 

0  0  0  -1 

10  0  0 
0-100 

and 

10  0  0 

0  10  0 

0  0-10 
0  0  0  -1 
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The  Dirac  equation  is  then  a  set  of  four,  coupled,  first  order  partial  differential 
equations  in  space  and  time.  The  four-component  wave  function  solution  to  this  equation 
corresponds  to  two  positive  energy  components  and  two  negative  energy  components, 
each  with  a  spin-up  and  spin-down  component  (Messiah,  1999:  888-892,920-924).  The 
a,  operators  are  velocity  operators,  while  fi  is  an  parity  operator.  Observable  quantities 

associated  with  this  new  internal  degree  of  freedom  are  t  he  energy,  relativistic  mass, 
current  density,  total  angular  momentum,  spin,  and  parity.  The  operators  associated  with 
these  observables  are  (Messiah,  1999:  921) 

energy:  H  =  e<j>  +  a-(p - )  +  j3m0c2, 

c 

(19) 

where  H  is  the  Dirac  Hamiltonian,  (f>  and  A  are  external  electric  and  magentic  potentials, 
respectively.  The  relativistic  mass  is  given  by 

M  =  H-e</>  (20) 

Current  density  is  found  using 

j(r)  =  a-S{f  -r0),  (21) 

where  S{f  -  r0  )  is  the  Dirac  delta  function.  Electron  spin  is  given  by 

S,=\&t  (22) 

where  er  is  the  ith  Pauli  spin  matrix.  Finally,  the  parity  of  the  wave  function  is  given  by 

P  =  J3P0,  (23) 
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where  Po  is  the  initial  parity.  These  Pauli  spin  matrices  can  be  expressed  in  terms  of  the 


aj  matrices  (Messiah,  1999:  891): 


and 


5,  =  -iaxay, 
5 x  =  -iaya, , 


oy=-iazax. 


(24) 

(25) 


(26) 


In  the  presence  of  an  external  field,  the  Dirac  Hamiltonian,  Hd,  becomes 

£>A  „ 

H D  =  e(f>  +  ca  -(p - )  +  j3m0c2.  (27) 

c 

For  the  hydrogen  atom,  in  the  absence  of  an  external  magnetic  field,  this  equation 
reduces  to 


fPP  =  [ c(a  ■  p )  +  j3m0c2  +  ecjPfi!  . 


(28) 


While  it  is  possible  to  construct  an  exact  solution  to  this  equation  in  terms  of 
spherical  harmonics  for  the  angular  coordinates  and  hypergeometric  functions  for  the 
radial  coordinate,  such  a  construction  does  not  shed  much  light  on  the  nature  of  the 
bound  energy  states.  The  details  of  the  solution  can  be  found  in  various  sources 
(Messiah,  1999:  930-933;  Balasubramanian,  1997:  159-175;  Bethe  et  al,  1957:  63-71). 
The  electronic  energy  levels  for  the  Dirac  hydrogen  atom  are  given  by  (Bethe  et  al,  1956: 
67-68) 


nj 

r  \ 

i+ 

Za 

i 

\n~j  +  \  +  'l{j+  iY  ~Z2a2  , 

(29) 
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The  first  term  is  the  non-relativistic  energy  for  the  bound  electronic  states  of  the 
hydrogen  atom.  Higher  order  corrections  involve  both  the  principle  quantum  number  n, 
as  well  as  the  total  angular  momentum  quantum  number  j.  This  illustrates  the  importance 
of  a  relativistic  picture  of  the  atom.  Corrections  to  the  non-relativistic  energy  increase 
roughly  as  Z4 .  Note  that  this  Taylor  series  expansion  in  powers  of  (Zaf  is  appropriate 
for  Za  «  1 .  This  expansion  leads  to  the  Russell-Sanders  spin-orbit  coupling  scheme. 
Such  an  approximation  is  not  valid  for  uranium,  where  Za  -  0.6 .  In  this  case, 

—  «  [Za  f  ,  and  the  electrostatic  electron-electron  interaction  can  be  treated  as  a 


perturbation  to  the  magnetic  interaction  between  the  electron  and  the  field  of  the  nucleus. 
This  approximation  leads  to  the  j-j  spin-orbit  coupling  scheme,  which  is  more  appropriate 
for  very  heavy  elements. 
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Detailed  examination  of  the  negative  energy  component  solutions  to  the  Dirac 
equation  for  the  free  electron  shows  in  the  non-relativistic  limit  where  E  -  m0c 2  «  m0c 2 , 
are  much  larger  than  the  negative  energy  components,  especially  in  the  valence  region 
(Balasubramanian,  1997:  143-144).  Thus,  the  four-component  Dirac  wave  function 
naturally  separates  into  two  large  and  two  small  components.  Rewriting  the  Dirac 
equation  in  terms  of  two,  coupled  differential  equations  with  two,  two-component  wave 
functions  yields  the  Pauli  approximation  to  the  Dirac  Hamiltonian  in  the  absence  of  an 
external  magnetic  field  (Balasubramanian,  1997:  145-147) 


=  E  +  e<p  +  -1- V2  +  +  elpf  +i-^-E-p 

2mn 


2  m0c 


2  m0c 


Mo 


2  m0c 


a  -(e  x  p)-  /J0(a  •  H ) 


(33) 


where  H  is  the  magnetic  field,  E  is  the  electric  field,  and  /u0  is  the  Bohr  magneton, 


defined  by 


Mo  = 


efi 

2  m0c 


(34) 


While  somewhat  cumbersome,  the  separate  terms  have  simple  interpretations. 
The  first  three  terms  are  the  non-relativistic  Schrodinger  Hamiltonian.  The  next  term  is 
the  mass-velocity  correction  that  accounts  for  the  variation  in  electron  mass  with  speed. 
The  fifth  tenn  is  known  as  the  Darwin  term,  and  is  a  result  of  “zitterbewegung”,  or 
trembling  motion.  It  is  a  result  of  the  Heisenberg  uncertainty  principle.  Non- 
relativistically,  the  uncertainty  in  the  location  of  an  electron  can  be  measured  to  any 
accuracy  using  higher  and  higher  energy  photons.  Relativistically,  there  is  a  limit  to  this 
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photon  energy  used  to  locate  the  electron,  because  at  photon  energies  above  2 m0c2 ,  pair 

production  can  occur.  This  results  in  an  effective  smearing  of  the  charge  of  the  electron 
(Balasubramanian,  1997:  186).  The  final  two  terms  account  for  the  spin-orbit  coupling 
between  the  intrinsic  electron  magnetic  moment  and  the  orbital  angular  momentum.  The 
successes  of  the  Dirac  equation  is  the  prediction  of  electron  spin  as  an  observable 
property  in  the  non-relativistic  limit,  as  well  as  accounting  for  the  correct  value  for  the 
electron  magnetic  moment.  Thus,  the  inclusion  of  electronic  spin  in  the  non-relativistic 
theory  as  an  additional  assumption  is  validated  and  explained  in  the  non-relativistic  limit 
of  the  Dirac  equation. 

Relativistic  Many-Electron  Hamiltonians 

Now  that  a  Lorentz-covariant  electronic  wave  function  is  available  for  the 
hydrogen  atom,  the  next  logical  step  is  to  try  to  extend  this  approach  to  larger  atoms  and 
molecules. 

A  relativistic  wave  function  for  a  many-electron  atom  can  be  constructed  as  the 
sum  of  the  one-electron  Dirac  Hamiltonians  along  with  an  electron-electron  interaction 
term 

H  =  'LK+'L‘>II,  (35) 

i  i<j 

where  the  ith  one-electron  Dirac  Hamiltonian,  hld  is  given  by 

hld  =  e<j)  +  ca  ■{pj - )  +  /3m0c~ ,  (36) 

c 

while  By  represents  a  general  electron-electron  interaction  term. 
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This  approach  is  analogous  with  the  non-relativistic  approximation  to  the  many-electron 
Hamiltonian  based  upon  the  one-electron  Schrodinger  equation  for  the  hydrogen  atom. 

Next,  a  Lorentz-covariant  description  of  electron-electron  interactions,  By,  is 
required.  Unfortunately,  the  electron-electron  interaction  requires  the  more  detailed 
treatment  afforded  by  quantum  electrodynamics,  where  vacuum  interactions,  virtual 
photon  exchanges,  and  electron  self-interactions  are  treated  perturbatively.  Even  then, 
there  is  no  closed  form  for  a  Lorentz-covariant  electron-electron  interaction.  Such 
interactions  must  be  treated  approximately  (Balasubramanian,  1997:  180;  Messiah, 
1999:  955-956;  Bethe  et  al,  1956:  170).  The  first  such  approximation,  widely  used,  is 
the  approximation  that  relativistic  corrections  to  the  electron-electron  interaction  are 
small  and  negligible,  and  that  the  Coulomb  interaction  is  an  appropriate  description, 
correct  to  zeroth  order.  This  leads  to  the  Dirac-Coulomb  Hamiltonian  (Kellogg,  1997: 
15) 


i<j  4^0  rij 


(37) 


Here,  r  is  the  interelectron  distance.  This  Hamiltonian  is  not  Lorentz-covariant; 


however,  corrections  to  the  Columbic  interaction  are  small  for  large  electron  separation, 
and  the  Dirac-Coulomb  Hamiltonian  is  quite  successful.  Another  approach  to 
determining  the  electron-electron  interaction  is  to  perturbatively  expand  the  quantum 
electrodynamics  interaction  term  in  powers  of  the  fine  structure  constant,  and  retain  those 
tenns  of  order  a 2 .  This  yields  the  Dirac-Coulomb-Breit  Hamiltonian  (Bethe  et  al,  1956: 
170;  Balasubramanian,  1997:  180;  Jackson,  1975:  593-595;  Breit,  1932:  616-624) 
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This  approximate  Hamiltonian  is  still  not  Lorentz-covariant;  however,  it  accounts 
for  most  of  the  chemically  relevant  electron-electron  interaction  effects.  The  first  tenn  is 
the  electrostatic  Coulomb  interaction  between  two  electrons,  the  second  term  accounts 
for  first  order  magnetic  interactions  between  the  intrinsic  magnetic  moments  of  the 
electron.  The  last  term  accounts  for  the  retardation  of  the  propagation  of  the 
electromagnetic  field  of  the  electron  due  to  the  finite  speed  of  light. 

Another  method,  based  on  perturbative  expansion  of  the  Dirac-Coulomb-Breit 
Hamiltonian  in  powers  of  of  yields  the  Breit-Pauli  Hamiltonian,  given  by 
(Balasubramanian,  1997:  193-194) 
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Here,  //  is  the  electron  magnetic  moment,  and  ?.  is  the  spin  angular  momentum 

of  the  ith  electron.  The  first  tenn,  Hnr,  is  the  non-relativistic  many-electron  Hamiltonian. 
The  second  tenn,  Hmv,  is  the  mass-velocity  term,  which  corrects  for  the  relativistic 
variation  in  electron  mass  near  the  speed  of  light.  The  third  term,  Hretardation,  corrects  for 
the  finite  propagation  speed  of  the  electromagnetic  field  of  the  electron.  The  fourth  term, 
HDarwin,  is  the  Darwin  correction,  described  earlier.  These  four  terms  comprise  scalar 
relativistic  effects,  and  do  not  require  a  two-component  wave  function  to  implement.  The 
fifth  term,  Hso,  is  the  spin-orbit  coupling  between  the  intrinsic  spin  angular  momentum 
of  the  electron  and  the  orbital  angular  momenta  of  all  the  electrons.  The  next  tenn,  Hss, 
is  the  spin-spin  coupling  between  the  intrinsic  spin-angular  momenta  of  multiple 
electrons.  In  order  to  incorporate  these  terms,  a  two-component  wave  function  is 
required.  The  last  tenn,  Hextemai,  involves  the  interaction  with  an  external  electric  and 
magnetic  field.  While  only  a  perturbative  treatment,  valid  for  light  atoms  and  molecules, 
the  Breit-Pauli  Hamiltonian  sheds  some  light  on  the  expected  effects  present  in 
relativistic  many-electron  Hamiltonians. 

With  a  well  defined,  albeit  approximate,  many-electron  Hamiltonian,  the  next  step 
in  constructing  relativistic  wave  functions  is  based  on  the  Hartree-Fock  mean  field 
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theory.  This  provides  the  first  theoretical  method  to  many-electron  systems;  however, 
electron  correlation  is  not  explicitly  included. 


Dirac  Hartree-Fock  Theory 

The  Hartree-Fock  (HF)  self-consistent  field  (SCF)  method  provides  the 
theoretical  framework  to  determine,  non-relativistically,  the  properly  antisymmetrized 
many-electron  single-detenninant  wave  function  for  atoms  and  molecules  (Szabo  et  al, 
1989:  108-152;  Levine,  2000:  305-312).  It  also  provides  the  basis  for  correlation 
calculations  through  multi-configuration  and  perturbation  methods. 

The  non-relativistic  many-electron  Hamiltonian  is  given  by  (Levine,  2000:  305) 

fl2  7n2  "  1  "  p1 

h—j- <47> 
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The  goal  of  the  HF  approximation  is  to  find  a  set  of  spin-orbitals,  which  minimize  the 
ground  state  electronic  energy.  It  is  a  variational  theory,  in  that  the  exact  ground  state 
energy  of  an  atom  or  molecule  is  a  lower  bound  to  the  HF  energy.  Additionally,  HF 
theory  is  a  single  determinant  theory.  This  means  that  the  ground  state  wave  function 
obtained  from  the  variationally  optimized  set  of  spin  orbitals  contains  only  a  single 
electron  configuration.  The  optimized  set  of  spin-orbitals  satisfy  the  equations  (Szabo  et 


al,  1989:  111-112) 
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where  %a  is  the  spatial  wave  function  of  the  ath  electron.  The  total  wave  function  is  the 


product  of  individual  electron  spatial  functions, 
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This  final  tenn  results  in  two,  one-electron  operators,  J  and  K.  The  first  operator 
is  the  Coulomb  operator  and  represents  the  average  electric  field  due  to  electron  two 
experienced  by  electron  one.  This  is  the  origin  of  the  mean-field  concept.  Each  electron 
experiences  an  average  potential  due  to  all  the  other  electrons.  The  second  term  is  an 
exchange  potential  arising  from  the  Pauli  Exclusion  Principle.  Because  the  total  wave 
function  must  be  antisymmetric  with  respect  to  the  exchange  of  two  electrons,  the  motion 
between  electrons  with  parallel  spins  is  correlated  in  the  HF  theory.  As  a  result,  electrons 
experience  an  exchange  potential,  quantum  mechanical  in  nature,  which  repels  electrons 
with  parallel  spins  and  prevents  them  from  occupying  the  same  orbital  (Szabo  et  al,  1989: 
111-115).  In  its  eigenvalue  form,  the  HF  equations  are 
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where  sa  is  the  energy  of  the  ath  electron.  The  one-electron  operators  are  defined  by 
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The  term  in  parentheses  of  the  first  equation  is  the  Fock  operator,  / (q  )  .  In 
operator  form,  this  equation  is 


f(rl)w,(r1)  =  £aWi(r1)- 


(56) 


Expanding  the  unknown  molecular  orbitals,  y/i  (r  )  ,  in  terms  of  a  finite  set  of 
known  basis  functions,  <j)k ,  yields 

Mfi)  =  Xc,,4  .  (57) 

i=l 

Here,  the  cu  are  molecular  orbital  coefficients.  This  reduces  the  HF  integro- 

differential  equation  to  a  set  of  algebraic  equations:  the  Roothan  equations  (Szabo  et  al, 
1989:  136-138) 

Ev^.Zv^  (58) 


where  F  is  the  Fock  matrix,  defined  by 
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(59) 


!)/(!¥,  (1), 

and  S  is  the  overlap  matrix,  defined  by 

Srr=\d?£m.(\).  (60) 

Thus,  finding  the  optimal  orbitals  that  minimize  the  Hartree-Fock  energy  consists 
of  solving  the  Roothan  equations  in  a  self-consistent  manner. 

The  non-relativistic  HF  theory  and  relativistic  Dirac-Hartree-Fock  (DHF)  theory 
are  analogous.  In  DHF  theory  (Saue  et  al,  1997:  937-948;  Oreg  ,  1975:  830-841; 
Aoyama  et  al,  1980:  1329-1332;  Matsuoka  et  al,  1980:  1320-1328;  Kim,  1967:  154-39; 
Lee  et  al:  1977:  5861-5876;  Dyall  et  al,  1991:  2583-2585),  the  non-relativistic 
Hamiltonian  is  replaced  with  the  Dirac-Coulomb  Hamiltonian,  and  the  spin-orbitals  have 
four  components  instead  of  one.  These  spinors  can  be  complex,  unlike  the  non- 
relativistic  case,  where  the  spin-orbitals  were  real.  This  four-component  wave  function  is 
also  expanded  in  a  real  basis,  as  was  done  in  the  non-relativistic  case,  where  (Saue,  1997: 
939) 
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Here,  L  and  S  represent  the  large  and  small  components,  respectively.  There  are  two 
sets  of  expansion  coefficients  now,  c“  and  cf  for  the  spin-up  and  spin-down 
components,  respectively.  In  a  manner  analogous  to  the  non-relativistic  theory,  these 
four-spinors  are  varied  until  the  total  electronic  energy  is  minimized.  This  leads  to  a 
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matrix  equation  similar  to  the  Hartree-Fock  theory,  except  that  the  Fock  matrix  and 


expansion  coefficients  are  now  complex.  The  Fock  matrix  splits  into  two  parts,  a  one- 
electron  matrix  (Saue,  1997:  940) 
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<67> 

As  in  the  non-relativistic  case,  V  represents  the  potential  energy  of  the  electron.  Atomic 
units  are  used  here,  where  h  =  e  =  m()  =  1 ,  in  order  to  simplify  the  expressions. 


The  two-electron  Fock  matrix,  is  given  by  (Saue,  1997:  940): 
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Here,  JLa  ,JLp  represents  the  Coulomb  operator  of  the  spin-up  ( a  )  and  spin-down  (/? ) 


electrons  respectively,  while  JSa  ,Jsp  represents  the  mean-field  due  to  positrons  with 
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spin-up  and  spin-down,  respectively.  The  exchange  operator,  K,  becomes  more 
complicated,  reflecting  the  possibility  of  exchange  between  electronic  ( L  )  components 
and  positronic  ( .S  )  components  with  various  spins. 

One  important  difference  exists  between  the  non-relativistic  Hartree-Fock  and  its 
relativistic  counterpart.  In  the  non-relativistic  theory,  the  energy  eigenvalue 
corresponding  to  the  Fock  operator  was  guaranteed  to  be  greater  than  or  equal  to  the 
exact  energy  by  the  variational  principle.  In  the  relativistic  case,  the  existence  of 
positronic  negative  energy  solutions  means  that  the  DHF  energy  is  not  bounded  from 
below.  In  the  non-relativistic  case,  there  were  no  negative  energy  solutions.  In  the 
relativistic  case,  electronic  solutions  look  like  excited  positronic  states,  and  unless  care  is 
taken  during  the  solution  of  the  relativistic  Roothan  equations,  variational  collapse  can 
occur.  This  occurs  because  of  the  fact  that  a  bound  electronic-positronic  state  is 
degenerate  with  an  unbound  electronic-postronic  state  in  the  relativistic  theory.  Thus, 
instead  of  variationally  optimizing  the  orbitals  as  in  the  non-relativistic  case,  the  orbitals 
are  minimized  with  respect  to  electronic  states  and  simultaneously  maximized  with 
respect  to  positronic  states.  Another  complication  with  the  relativistic  DHF  theory  is  that 
the  basis  sets  for  the  large  and  small  components  are  related  via  a  “kinetic  balance” 
requirement  (Dyall,  1991:  2585) 

Xs  =  ^  ■  VzL  •  (69) 

Here,  the  small  component  (positronic)  wave  functions  are  related  to  the  2x2 
Pauli  spin-matrix,  <r ,  operating  on  the  gradient  of  the  large  component  (electronic)  wave 
function. 
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DHF  methods  suffer  several  computational  difficulties  (Saue,  1997:  937-938). 
First,  spin  and  spatial  wave  functions  are  coupled,  resulting  in  complex  total  wave 
functions  in  general.  Thus,  spatial  and  spin  symmetry  can  not  be  handled  separately, 
which  has  resulted  in  substantial  computational  savings  for  non-relativistic  computations. 
Another  computational  difficulty  lies  in  the  fact  the  DHF  basis  sizes  are  generally  much 
larger  than  their  non-relativistic  counterparts.  Both  the  large  and  small  components  are 
expanded  in  separate,  but  coupled  basis  sets.  The  small  basis  set  size  can  be  generally 
twice  the  size  of  the  large  component  basis.  Computationally,  the  Hartree-Fock  method 
scales  minimally  as  N  ,  where  N  is  the  number  of  basis  functions.  Thus,  Dirac-Hartree- 
Fock  calculations  typically  involve  an  order  of  magnitude  or  more  increase  in 
computational  complexity  over  non-relativistic  Hartree-Fock  computations.  These 
difficulties  currently  limit  DHF  methods  to  atoms  and  some  small  molecules.  A  recent 
DHF  calculation  (de  Jong,  et  al,  1999:  45)  for  U022+ compares  the  non-relativistic  and 
relativistic  results  for  the  ground  electronic  states.  Electron  correlation  was  included  in 
this  calculation  via  coupled  cluster  singles  and  doubles  with  some  triples  (CCSD(T)) 
(Raghavachari,  1996:  12964-12965),  and  this  calculation  represents  perhaps  the  all¬ 
electron  computational  state-of-the-art  on  the  ground  state  of  uranyl  (de  Jong,  1999:  41- 
52).  Figure  3  reproduces  the  electronic  states  from  the  paper. 
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Figure  3.  Electronic  Ground  State  for  Uranyl  Using  Both  Non-relativistic  HF  and  Four- 
component,  Fully  Relativistic  DFIF  Methods  (de  Jong,  1999:  45) 

Dirac-Hartree-Fock  theory  can  result  in  fairly  accurate,  fully-relativistic 
calculations  of  molecular  electronic  ground  states.  Like  the  non-relativistic  counterpart, 
DHF  provides  the  best,  single-determinant  wave  function,  ignoring  electronic  correlation 
effects.  As  such,  HF  and  DHF  can  describe  only  the  electronic  ground  state.  In  order  to 
describe  the  excited  electronic  states,  electronic  correlation  needs  to  be  incorporated,  and 
the  wave  function  must  be  expanded  in  a  series  of  determinants.  As  was  evident  in  the 
DHF  equations,  the  fully  relativistic  treatment  can  become  quite  complicated  and 
computationally  intensive.  What  is  needed  is  a  method  that  is  a  compromise  between  the 
non-relativistic  and  fully-relativistic  methods  that  also  provides  computational  savings 
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when  computing  electron  correlation  energies  using  multi-configuration  methods. 
Relativistic  effective  core  potentials  prove  to  be  just  such  a  compromise. 

Relativistic  Effective  Core  Potentials 

The  DHF  method  provides  a  basis  for  constructing  the  approximate  relativistic 
wave  function  for  atoms  and  molecules.  It  is  an  all-electron  method,  and  the  number  of 
electrons  that  must  be  treated  explicitly  for  actinide  elements  and  molecules  is  large. 
Given  the  requirements  for  kinetic  balance  between  the  basis  sets  for  the  large  and  small 
components,  the  number  of  basis  functions  required  for  accurate  treatment  of  heavy 
atoms  and  molecules  can  easily  become  formidable  for  interesting  actinide  molecules. 
Thus,  a  method  for  accurately  treating  relativity  and  only  explicitly  including  those 
electrons  that  are  relevant  chemically  is  attractive.  The  relativistic  effective  core 
potential  (RECP)  method  accomplishes  both  of  these  goals. 

Several  different  types  of  RECPs  exist.  Among  them  include  potentials  by 
Kuchle,  Dolg,  Stoll,  and  Preuss  (Kiichle  et  al,  1994:  7535-7542),  Hay  and  Wadt  (Hay, 
1993:  5469),  and  Christiansen  style  RECPs  (Ennler  et  al,  1988:  139-182),  among 
others.  The  Christiansen  style  RECPs  were  selected  for  this  work  for  their  spin-orbit 
operator. 

The  RECP  method  (Pacios,  1985:  2664-2671;  Lee  et  al,  1977:  5861-5876; 
Ermler  et  al,  1991 :  829-846;  Ennler  et  al,  1988:  139-182)  is  based  upon  the  assumption 
that  atomic  or  molecular  electrons  can  be  separated  into  an  inert  core  region  and  a 
chemically-active  valence  region.  The  RECPs  are  generated  from  the  large  component 
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valence  electron  solutions  to  numerical  DHF  calculations  of  the  form  (Pacios  et  al,  1985: 


2664) 


f  -  v2 


-  —  +  U™CP(r)  +  Jli-Kli 


xij  =  £ijX,j  ■ 


(70) 


Here,  Zv  is  the  number  of  valence  electrons,  while  J and  Ktj  represent  the 

Coulomb  and  exchange  tenns.  Atomic  units  are  used  here.  This  equation  is  inverted  in 
order  to  recover  the  core  potential  term  for  a  particular  value  of  l  and  j  .  The  radial 

dependence  of  the  pseudospinor,  Xu ,  is  determined  from 


xii(r)  =  ViM)  +  FiM) 


(71) 


Xy  is  referred  to  as  a  pseudospinor,  because  it  approximates  the  two-component, 

electronic  portion  of  the  fully-relativistic  four-component  spinor. 

The  wave  functions,  y/t. ,  are  determined  from  the  large  component  DHF  solutions, 

while  the  second  term,  FJr)  is  selected  to  cancel  radial  oscillations  in  the  core  region 

and  eliminate  nodes  in  the  pseudospinor.  It  is  critical  to  generate  a  nodeless 
pseudospinor  in  order  to  invert  equation  70  and  recover  the  core  potential  (Pacios,  1985: 
2664).  The  core  potentials  generated  are  of  the  fonn  (Ermler  et  al,  1991 :  829-830) 


jjRECP  _  jjRECP 


L  1  M  j 

0)  +  X  Z  - URLjCP{r •)) •  Ijrrij }(/ym/. |  (72) 


1=0  j=i-\  mj=-j 


Here,  I  Ijm^j  and  (ijnij  |  are  projection  operators,  ensuring  the  core  potentials 

interact  with  those  parts  of  the  wave  function  with  the  corresponding  orbital  angular 
momentum,  / ,  total  angular  momentum,  j ,  and  z-component  of  total  angular 
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momentum,  .  While  the  sum  over  /  is  infinite  in  principle,  truncation  of  the  sum  at  an 

orbital  angular  momentum  one  greater  than  the  largest  /  quantum  number  of  the  core 
electrons  introduces  negligible  error  (Ermler  et  al,  1991 :  830).  A  weighted  average 
relativistic  effective  core  potential  is  then  generated  by  (Ermler  et  al,  1991 :  830;  Ermler 
etal,  1999:  152-153) 

UARECP = Uf  RECP(r)  +  £  YjyUARECP(r)  -  UfECP(r))  ■  \  Im,  )(lm,  | ,  (73) 
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and 
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The  purpose  for  generating  the  average  relativistic  effective  core  potential  is  to  separate 
the  scalar  relativistic  effects  (mass-velocity,  Darwin  tenns)  from  the  spin-orbit  effects. 
Scalar  relativistic  effects  can  be  incorporated  into  spin-free  wave  functions.  Spin-orbit 
effects  require  the  introduction  of  electron  spin,  and  consequently,  require  a  two- 
component  wave  function.  An  effective  spin-orbit  Hamiltonian  is  generated  by  (Ermler 
etal,  1991:  830;  Ermler  et  al,  1988:  152) 

=  S  ■j^-^—AUf‘cr(r)Y,\lml)(lmlf.\lmr){lmi\,  (75) 
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where 

A U™CP  (r)  =  U™ c{  (r)  -  (/  + 1)  •  C/REfE  (r) .  (76) 

This  spin-orbit  Hamiltonian  can  be  written  in  the  form  (Ermler  et  al,  1988:  153) 

\Im,)(lm,\.  (77) 
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The  function  g{r )  is  a  constant  that  depends  on  the  radial  quantum  number  of  the  orbitals 


involved.  The  projection  operators  \lml'j(lml  |  in  all  forms  of  the  effective  core  potentials 

ensure  that  orbitals  with  the  correct  angular  momentum  interact  with  the  correct  core 
potential  tenn.  The  spin-orbit  potential  fonned  this  wave  acts  as  a  one-electron  operator. 
As  such,  it  is  only  approximately  valid  for  the  Dirac-Coulomb  Hamiltonian.  Higher- 
order  electron-electron  interactions  represented  by  the  Dirac-Coulomb-Breit  Hamiltonian 
are  not  effectively  modeled  using  this  spin-orbit  operator. 

Once  the  core  potential  and  spin-orbit  potentials  are  obtained,  they  are 
approximated  by  Gaussian  function  fits  of  the  fonn 

1  M 

UfRECP{r )  =  -VC/%eV  ,  (78) 

r  ,= i 

with  expansion  coefficients  Cu  and  exponents  «/; ,  facilitating  their  use  in  standard  ab 
initio  computational  packages. 

The  benefit  of  RECPs  is  that  they  can  be  used  to  incorporate  scalar  relativistic 
effects  in  normally  non-relativistic  theoretical  methods  such  as  HF,  MCSCF,  or  Density 
Functional  Theory  (DFT).  In  addition,  RECPs  offer  a  way  to  substantially  reduce  the 
number  of  explicitly  treated  electrons.  RECPs  are  currently  the  most  commonly  used 
method  for  theoretical  studies  of  heavy  element  molecules. 

Electron  Correlation  Models 

The  development  of  accurate  relativistic  core  and  spin-orbit  potentials  enables  the 
inclusion  of  relativistic  effects  in  non-relativistic  theoretical  methods.  This  allows  the 
successful  application  of  powerful  theoretical  methods  developed  for  use  in  light 
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molecule  calculations  to  be  applied  to  heavy  element  molecules  at  a  reasonable 
computational  cost.  While  this  section  will  discuss  electron  correlation  techniques  from 
a  non-relativistic  standpoint,  all  have  either  been  applied  to  relativistic  calculations 
directly,  as  with  MP2  and  MCSCF,  or  indirectly,  via  RECPs  incorporated  into  the  non- 
relativistic  technique. 

Electronic  correlation  effects  are  not  directly  observable,  according  to 
Raghavachari  (Raghavachari,  1996:  12960).  Electronic  correlation  is  a  measure  of  that 
energy  that  is  not  accurately  modeled  by  the  Hartree-Fock  theory.  In  fact,  correlation 
energy  is  defined  as  the  difference  between  the  exact  non-relativistic  energy  and  the 
Hartree-Fock  energy.  This  inaccuracy,  inherent  in  the  Hartree-Fock  model,  incorrectly 
models  the  dissociation  of  molecules  into  constituent  atoms.  Improvements  to  HF  that 
model  this  energy  accurately  can  predict  the  correct  dissociation  energies  for  molecules, 
among  other  properties. 

Since  the  HF  single  determinant  wave  function  frequently  accounts  for  a  large 
fraction  of  the  total  electronic  energy,  one  popular  and  efficient  technique  for 
determining  electron  correlation  energies  is  through  perturbation  theory.  The  electron 
correlation  is  treated  as  a  small  perturbation  to  the  HF  wave  function,  and  corrections  to 
the  HF  wave  function  and  energy  are  computed  via  a  perturbation  expansion.  The  most 
frequently  used  is  second  order  perturbation  theory,  also  known  as  Moller-Plesset  second 
order  perturbation  theory  (MP2).  In  MP2,  the  Fock  operator  is  the  zeroth-order 
Hamiltonian,  and  both  the  HF  wave  function  and  HF  energy  are  expanded  in  a  power 
series  of  the  perturbation  (Raghavachari  et  al,  1996:  12962;  Szabo  et  al,  1989:  350-353). 
Using  this  perturbative  technique,  one  can  show  that  the  HF  energy  is  correct  to  first 
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order,  and  corrections  to  this  energy  occur  at  second  order  and  beyond.  Higher  order 
perturbation  corrections  are  possible;  however,  MP2  is  the  most  frequently  used.  The 
perturbed  Hamiltonian,  H' ,  is  given  by 

H'=H0+V.  (79) 

The  zeroth-order  Hamiltonian  is  the  Fock  operator, 

H>='Zki  +  'ZJt-K9’  <80> 

i  i<j 

while  the  perturbation  is  given  by 

v  =  'Lt-ji  +  ki-  <81> 

i<j  ij 


The  sum  of  the  zeroth  and  first-order  energies  correspond  to  the  HF  energy.  The  second 
order  correction  to  the  energy  is  given  by  (Szabo  et  al,  1989:  351). 


^2>  =  £ 


Xi  V  Xn 


E(0)  -E 


(0) 


(82) 


Here,  the  unperturbed  N  electron  wave  function  is  written  as  the  product  of  N  spatial 
wave  functions, 

klT  —  y  y  v  y  .  .  .  y 

T0  —  /C0A1/C2  •••An  AN  (§3) 

and  E'f  '  represents  the  Hartree-Fock  energy  of  the  ith  spatial  function.  The  result  of  this 

second-order  perturbation  expansion  of  the  Hartree-Fock  wave  function  is  to  introduces 
some  single  and  double  excited  state  energies  into  the  ground  state  energy.  Singly 
excited  wave  functions  do  not  couple  directly  with  the  HF  wave  function.  They  do; 
however,  couple  to  the  doubly  excited  wave  functions,  which  in  turn  are  coupled  to  the 
HF  wave  function.  This  is  proven  in  Brillouin’s  Theorem  (Szabo  et  al,  1989:  128-129). 
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MP2  theory  scales  computationally  as  the  number  of  basis  functions  (K)  to  the 
fifth  power  (K5)  (Raghavachari  et  al,  1996:  12962),  making  it  attractive  compared  to  HF 
alone,  which  scales  as  K4.  Higher  order  perturbations  scale  as  the  sixth,  seventh,  and 
eighth  power  for  third,  fourth,  and  fifth  order  perturbations  respectively,  making 
perturbations  beyond  fourth  order  computationally  demanding  and  seldom  used. 

Another  common  technique  for  modeling  electron  correlation  energy  is  through  a 
multi-configurational  approach,  which  involves  expanding  the  wave  function  in  a  linear 
combination  of  excited  electronic  configurations  (Raghavachari  et  al,  1996:  12962). 
Various  expansion  terms  then  represent  exciting  electrons  from  the  reference 
configuration  into  a  number  of  excited  configurations.  Then  all  that  is  necessary  is  to 
determine  the  expansion  coefficients.  If  all  possible  configuration  state  functions  were 
used,  the  expansion  would  be  exact;  however,  in  practice,  this  expansion  must  be 
truncated.  As  with  the  finite  basis  set  expansion  for  the  single  determinant  wave 
function,  this  truncation  introduces  error  into  the  correlated  wave  function.  This  multi- 
configurational  technique  is  necessary  to  accurately  describe  excited  electronic  states. 
The  reason  for  this  is  that  excited  electronic  states  with  the  same  term  symbols  mix  with 
the  ground  state.  The  Hartree-Fock  neglects  this  mixing  by  writing  the  wave  function  as 
a  single  Slater  determinant,  however,  the  actual  ground  state  wave  function  contains 
contributions  with  excited  states  with  the  same  term  symbol.  A  more  accurate  wave 
function  expansion  requires  multiple  Slater  detenninants,  where  each  determinant 
represents  a  particular  state  (ground  and  excited). 

One  method  for  improving  upon  the  HF  wave  function  is  to  include  a  relatively 
small  number  of  additional  configurations  (Raghavachari  et  al,  1996:  12967).  This 
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technique  is  the  multi-configuration  self-consistent  field  method  (MCSCF).  In  the 
MCSCF  approach,  the  wave  function  is  expanded  in  a  small  set  of  configuration  wave 
functions  (Szabo  et  al,  1989:  258). 

^mcsc^^o+^+.-.  +  c,^,  (84) 

where  the  c,  are  Cl  coefficients.  Then,  both  the  Cl  and  molecular  orbital  coefficients  are 
variationally  optimized.  The  MCSCF  method  is  identical  to  the  FIF  method  if  only  one 
expansion  wave  function  is  used.  Reference  configurations  can  be  chosen  to  include 
chemically  relevant  states,  or  states  that  may  be  nearly  degenerate  with  the  HF  ground 
state.  Flowever,  this  approach  can  introduce  bias  into  the  calculation  (Raghavachari  et  al, 
1996:  12967).  A  way  to  avoid  this  bias  is  to  identify  a  set  of  active  orbitals,  typically  the 
highest  occupied  molecular  orbitals  (FIOMO),  and  include  all  possible  excited 
configurations  arising  from  this  active  space.  This  is  a  complete  active  space  (CAS) 
approach,  ft  is  beneficial  to  use  a  CAS  approach  when  the  excited  states  of  a  molecule 
may  not  be  well-known;  however,  the  number  of  configurations  included  in  a  CAS 
MCSCF  can  quickly  become  enormous  if  the  active  space  is  large,  or  there  are  a  large 
number  of  virtual  orbitals. 

Another  benefit  to  the  MCSCF  method  is  that  it  provides  an  improved  wave 
function  for  use  in  the  more  general  configuration  interaction  (Cl)  approach.  The  Cl 
wave  function  is  written  as 

y  =c»*0  +EcW  +  I‘f  >F;  +  +....  (85) 

a  a,b  a,b,c 

Where  the  first  term  is  the  HF  wave  function,  the  second  tenn  is  a  sum  over  all  singly 
excited  states  from  the  HF  ground  state,  the  third  term  consists  of  a  sum  over  all  doubly 
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excitations,  etc.  Including  all  possible  determinants  is  a  full  Cl  expansion,  which 
provides  the  correlation  exactly;  however,  such  an  expansion  is  almost  never  possible 
except  for  very  small  systems  using  minimal  basis  sets.  In  practice,  this  Cl  expansion 
must  be  truncated  at  some  point.  It  can  be  shown  (Szabo  et  al,  1989:  128-131)  that  there 
is  no  direct  mixing  between  the  HF  ground  state  and  the  singly  excited  states.  However, 
the  HF  ground  state  does  mix  with  the  doubly  excited  states,  which  in  turn  couples  to  the 
singly  excited  states.  Thus,  the  first  logical  place  to  truncate  the  Cl  expansion  in  order  to 
obtain  an  improvement  over  the  HF  wave  function  is  to  include  single  and  double 
excitations  (CISD).  The  CISD  wave  function  is  given  by 

fa® = <vr0 + .  m 

a  a,b 


where  the  expansion  coefficients  are  variationally  optimized.  The  CISD  method,  while 
conceptually  simple,  is  very  demanding  computationally.  In  an  N-electron  system  using 
K  basis  functions,  the  number  of  all  possible  double  excitations  is  given  by  (Szabo  et  al, 
1989:  234) 


'nVik-n 

v  ^ 


(87) 


Thus,  the  CISD  method  scales  roughly  as  (2K)~  N2 .  Truncating  the  expansion  at  the 
quadruple  excitation  level  scales  roughly  as  (2K)4  N4 ,  while  sextuple  excitations  scales 
as  (2K)6  N 6 .  It  is  apparent  that  the  Cl  method  increases  in  computational  difficulty  very 
quickly.  While  the  full  Cl  does  scale  linearly  with  the  size  of  the  system,  the  truncated 
Cl  expansions  do  not. 
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By  allowing  single  and  double  excitations  from  a  multi-reference  wave  function, 
as  opposed  to  the  single  reference  wave  function,  some  triple  and  quadruple  excitations 
are  effectively  included  in  the  CISD  expansion.  References  are  those  electronic 
configurations  from  which  electrons  can  be  excited  out  of  in  the  Cl  expansion.  This  is 
known  as  multi-reference  CISD  (MR-CISD).  Single-  and  multi-reference  Cl  expansion 
methods  are  some  of  the  few  theoretical  methods  available  for  exploring  atomic  and 
molecular  excited  states. 

One  important  point  to  note  when  using  Cl  expansion  methods  is  that  the  choice 
of  orbitals  can  impact  the  accuracy  of  the  Cl  expansion.  There  are  two  types  of  orbitals 
one  can  use  when  performing  a  Cl  expansion:  molecular  and  natural  orbitals.  Natural 
orbitals  are  those  orbitals  which  diagonalizes  the  density  matrix,  /(x^x^) ,  given  by 
(Szabo  et  al,  1989:  252-257) 

r(x  i ,  *i ')  =  X  &  )  V;  <*/)’  (8g) 

V 

where  y  and  y,j  are  defined  by 

Yij  =  \  dx,  dx, '  x]  (i,  )y(xx ,  T, ' )Zj  (x, ' )  (89) 

The  Xi  in  equation  88  are  the  spatial  orbitals,  while  T'(xl  ,x2,---,xv  )  is  the  electronic 

wave  function,  which  is  the  product  of  the  spatial  orbitals. 

The  molecular  orbitals  diagonalize  the  Fock  matrix,  but  they  do  not  diagonalize 
the  density  matrix.  A  Cl  expansion  using  the  natural  orbitals  converges  faster  than  one 
formed  from  molecular  orbitals  (Szabo  et  al,  1989:  255).  Thus,  a  Cl  expansion  using 
molecular  orbitals  will  require  more  configuration  state  functions  to  achieve  the  accuracy 
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of  an  expansion  fonned  from  natural  orbitals.  In  many  cases,  the  difference  in  accuracy 
may  be  negligible;  however,  there  are  some  cases  were  this  accuracy  difference  can  be 
substantial.  One  difficulty;  however,  is  that  in  order  to  obtain  the  natural  orbitals,  one 
needs  the  Cl  wave  function  (Szabo  et  al,  1989:  257).  Luckily,  approximate  natural 
orbitals,  formed  from  a  small  number  of  configurations,  such  as  an  MCSCF  wave 
function,  are  almost  as  good  as  natural  orbitals  obtained  from  the  Cl  wave  function 
(Szabo  et  al,  1989:  257). 

One  final  note  for  electron  correlation  techniques  is  in  order,  especially  when 
incorporating  relativistic  spin-orbit  effects.  Because  electron  correlation  is  determined 
primarily  by  a  two-electron  operator  between  pairs  of  reference  configurations  and 
between  reference  configurations  and  doubly  excited  configurations,  it  is  generally  more 
difficult  to  converge  computationally  than  spin-orbit  effects.  Spin-orbit  effects  are 
expressed  in  the  core  potential  formulation  as  one-electron  operators  between  pairs  of 
reference  configurations  and  singly  excited  configurations  (Yabushita  et  al,  1999:  5792). 
As  a  result  of  the  two-electron  operator,  accurately  modeling  electron  correlation  is  a 
more  difficult  problem  than  inclusion  of  relativistic  spin-orbit  effects,  and  efficiently 
handling  the  electronic  correlation  will,  in  general,  lead  to  accurate  inclusion  of  spin- 
orbit  effects.  And  so,  when  incorporating  spin-orbit  effects  into  electronic  structure 
calculations,  it  makes  sense  to  include  them  in  conjunction  with  electron  correlation  at 
the  same  time,  rather  than  as  separate  steps  (Yabushita  et  al,  1999:  5792). 

While  the  previous  methods  for  modeling  electron  correlation  are  all  based  upon 
approximations  to  the  many-electron  wave  function  (Head-Gordon,  1996:  13218),  one 
method  which  is  not  based  on  wave  functions,  has  seen  tremendous  growth  and 
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application  in  the  past  decade.  This  technique  is  known  as  Density  Functional  Theory 
(DFT),  and  it  is  based  upon  the  premise  that  the  exact  molecular  ground  state  energy  is  a 
functional  of  the  electron  density  (Nagy,  1998:  5).  The  electronic  density,  p,  of  an  n- 
electron  system  is  a  function  of  three  variables,  while  the  n-electron  wave  function  is  a 
function  of  3n  variables.  This  reduction  in  degrees  of  freedom  results  in  substantial 
computation  savings  with  DFT  over  wave  function  based  methods.  This  energy  is  given 
by  (Head-Gordon,  1996:  13218) 

E(p)  =  T(p)  +  Ven  (p)  +  J(p)  +  Vxc  {p) .  (90) 

Fiere,  T  is  the  kinetic  energy,  J  is  the  Coulomb  repulsion,  Ven  is  the  electron- 
nuclear  interaction,  and  Vxc  is  the  electron  exchange  and  correlation  interaction.  The 
Coulombic  electron-electron  and  electron-nuclear  terms  are  straightforward  and  can  be 
computed  classically,  using  a  non-interacting  Fermi  gas  model  derivable  from  a  statistical 
treatment  of  the  atom  has  a  degenerate  electron  gas.  The  kinetic  energy  and  exchange- 
correlation  terms  cannot.  Kohn  and  Sham  showed  (Kohn  et  al,  1996:  12974-12975; 
Nagy,  1998:  7-12)  that  this  energy  can  be  recast  into  a  fonn  where  the  kinetic  energy  is 
approximated  by  a  system  of  non-interacting  electrons.  This  results  leads  to  the  Kohn- 
Sham  equations  (Nagy,  1998:  9) 

V2  'i 

-  ~  +  v(r )  +  Vj  (i r )  +  vxc  (?)  u,  =  £iui ,  (91) 

V  1  J 

Vj(?)  =  \^Xd?',  (92) 

j  \r  _r’| 

and 
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Here,  v(r )  is  an  external  potential  and  ut  is  the  wave  function.  Thus,  the  electron 

exchange  and  correlation  interactions,  as  well  as  the  fraction  of  kinetic  energy  not 
accounted  for  in  the  non-interacting  gas  model  are  wrapped  into  the  exchange-correlation 
functional.  Various  attempts  at  developing  an  exact  fonn  for  this  functional  have  been 
unsuccessful;  however,  several  approximate  functionals  have  been  developed  and 
successfully  applied  to  electronic  structure  calculations  (Head-Gordon,  1996:  13218). 
The  most  successful  correlation  functionals  have  been  gradient-corrected  approximate 
(GCA)  correlation  functionals,  while  successful  exchange  functionals  are  based  upon 
parameter  fits  to  the  exchange  energies  of  noble  gases  (Head-Gordon,  1996:  13218; 
Nagy,  1998:  41-45).  Combinations  of  the  GCA  and  semi-empirical  functionals  form 
hybrid  functionals,  which  also  have  been  very  successful  in  treating  many  ground-state 
molecular  properties.  The  accuracy  achieved  with  DFT  frequently  rivals  or  beats  MP2 
results,  with  less  computational  effort.  One  drawback  to  DFT  at  present;  however,  is  that 
DFT  is  a  ground  state  theory  only.  DFT  has  not  been  successful  at  describing  excited 
electronic  energies.  In  addition,  unlike  wave  function  based  theories  ,  there  is  no  known 
way  to  systematically  improve  the  results  obtained  with  DFT.  Without  the  ability  to 
accurately  describe  excited  electronic  states  directly,  or  approximate  excited  electronic 
state  energies  via  a  systematic  improvement  process  (perturbation  theory),  DFT,  despite 
its  successes,  is  simply  not  an  option  when  calculation  theoretical  electronic  spectra  of 


uranium  oxides. 


Summary 


This  chapter  briefly  covered  some  of  the  theory  behind  relativistic  electronic 
structure  calculations.  As  evidenced  by  results  from  perturbation  theory,  relativistic 
effects  must  be  included  for  heavy  element  compounds,  uranium  in  particular. 
Additionally,  electron  correlation  effects,  expected  to  be  large  in  actinide  molecules  pose 
additional  theoretical  challenges.  As  a  result,  theoretical  relativistic  electronic  structure 
calculations  on  actinide  molecules  have  pushed  the  state-of-the-art  in  ab  initio  molecular 
calculations  for  many  years,  and  they  will  continue  to  do  so. 

In  order  to  accurately  model  excited  electronic  states  in  uranium  oxides, 
relativistic  effects  must  be  incorporated,  and  some  type  of  multi-configurational  approach 
must  be  used.  Because  of  the  large  number  of  electrons  that  must  be  treated  in  order  to 
accurately  describe  excited  electronic  states,  a  fully-relativistic  method,  such  as  Dirac- 
Hartree-Fock,  coupled  with  either  many-body  perturbation  theory  or  a  multi- 
configuration  method  is  too  computationally  demanding.  By  using  a  relativistic  effective 
core  potential  method,  relativistic  effects  are  incorporated  to  first  order,  and  the  number 
of  electrons  that  must  be  treated  explicitly  in  advanced  correlation  techniques  is  reduced. 
These  compromises  prove  to  be  acceptable  in  the  chemically  active  valence  region  of 
uranium,  which  is  responsible  for  the  optical  fluorescence  spectra  in  uranyl,  as  well  as 
low  energy  photoelectron  spectra. 
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III.  Resources  and  Methodology 


The  complexity  involved  in  theoretical  studies  of  heavy-element  molecules 
requires  powerful  computing  platforms  and  specialized  software.  While  theoretical 
calculations  on  small,  light-element  molecules  can  be  performed  to  within  chemical 
accuracy  on  personal  computers  (PCs),  comparable  accuracy  with  actinide  molecules 
often  requires,  large,  massively  parallel  computers  with  large  amounts  of  memory,  disk 
storage  capacity,  and  high  input/output  throughput.  Also,  software  must  be  specifically 
designed,  and  theoretical  algorithms  carefully  constructed,  in  order  to  make  optimal  use 
of  these  computing  platforms. 

Hardware  Resources 

The  theoretical  calculations  performed  in  this  study  used  parallel  computing 
clusters  at  AFIT  and  at  the  Aeronautical  System  Center’s  (ASC)  Major  Shared  Resource 
Center  (MSRC). 

The  parallel  computing  cluster  at  AFIT  consists  of  64  PCs  running  Redhat  Linux 
7.3  operating  system.  Sixteen  of  these  PCs  have  1.4  GHz  Athlon  Advanced  Micro 
Devices  (AMD)  processors  with  775  megabytes  (Mb)  random  access  memory  (RAM),  20 
gigabyte  (Gb)  local  disk  space,  and  they  are  networked  via  fast  Ethernet.  The  other  48 
PCs  each  contain  dual  1  GHz  Pentium  III  processors,  1  Gb  RAM,  and  20  Gb  local  hard 
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disk  storage.  Thirty-two  of  these  dual-processors  PCs  are  networked  via  fast  Ethernet, 
while  the  other  16  dual-processor  PCs  are  networked  via  Myranet. 

Three  parallel  clusters  were  used  at  the  ASC  MSRC.  The  first  two  clusters 
consist  of  Compaq  ES  40  and  ES  45  machines.  The  ES  40  cluster  consists  of  16 
symmetric  multiprocessor  (SMP)  nodes,  with  each  node  containing  four,  833  MHz 
processors  and  4  Gb  shared  memory  per  node.  The  cluster  has  access  to  a  total 
workspace  of  about  800  Gb  (HP/Compaq  ES40/45  User’s  Guide,  2002:  4).  The  Compaq 
ES  45  cluster  consists  of  81  and  128  SMP  nodes  containing  four,  1  GHz  processors  per 
node.  Each  node  has  4  Gb  shared  memory,  and  the  cluster  has  a  total  workspace  of  8  Tb 
(HP/Compaq  ES40/45  User’s  Guide,  2002:  4). 

The  third  cluster  consisted  of  132  IBM  SP  3  SMP  nodes.  Each  node  contains 
four,  375  MHz  processors,  4  Gb  shared  memory,  and  2.4  terabytes  (Tb)  workspace  (IBM 
SP  P3  User’s  Guide,  2002:  1). 

Software  Resources 

This  research  project  used  two  software  packages.  The  first  package,  NWChem 
version  4.0.1  ,  was  developed  by  the  Molecular  Sciences  Software  group  of  the  Theory, 
Modeling  &  Simulation  program  of  the  Environmental  Molecular  Sciences  Laboratory 
(EMSL)  at  the  Pacific  Northwest  National  Laboratory  (PNNL).  It  was  designed  to 
perform  a  wide  range  of  molecular  calculations  on  massively  parallel  computing 
platforms  (NWChem,  2002;  High  Performance  Computational  Chemistry  Group,  2002). 
NWChem  was  used  for  geometry  optimizations,  vibrational  frequency  calculations,  and 

3  NWChem  home  page  is  at  http://www.emsl. pnl.gov:2Q8Q/docs/nwchem/nwchem.html 
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some  single  point  energy  calculations  without  spin-orbit  potentials  on  uranyl  and  uranate 
ions. 

NWChem  was  compiled  and  tested  on  the  AFIT  cluster  using  the  GNU  FORTRAN  77 
compiler.  It  was  also  installed  and  tested  on  the  MSRC  Compaq  ES40/45  clusters  by 
MSRC  personnel. 

The  second  software  package  used  was  COLUMBUS  version  5.8.1. 
COLUMBUS4  (Lischka,  2001;  Lischka,  2000;  Tilson,  2000;  Yabushita,  1999;  Sheppard, 
1988)  is  a  suite  of  LORTRAN  programs  developed  at  Ohio  State  University  and  Argonne 
National  Laboratories,  and  it  is  maintained  at  the  University  of  Vienna,  Austria. 
COLUMBUS  was  used  to  perform  SCF,  MCSCF,  and  two-component  MR-CISD 
calculations  on  both  the  uranyl  and  uranate  ions.  COLUMBUS  was  compiled  and  tested 
on  the  AFIT  cluster  using  the  GNU  FORTRAN  77  compiler.  MSRC  personnel  installed 
and  tested  COLUMBUS  5.8. 1  on  the  IBM  SP  P3  cluster.  The  SCF  and  MCSCF 
calculations  did  not  use  spin-orbit  potentials;  however,  the  MR-CISD  calculations 
included  spin-orbit  potentials  for  uranium  and  oxygen.  Unfortunately,  COLUMBUS 
5.8. 1  does  not  included  a  parallelized  version  of  the  spin-orbit  Cl  (SOCI)  program. 

While  this  had  no  effect  on  the  calculation,  it  limits  the  practical  number  of  configuration 
state  functions  (CSFs)  in  the  Cl  calculation  to  about  10  million.  COLUMBUS  5.9  does 
have  a  parallel  multi-reference  SOCI  (MR-SOCI)  program. 


4  COLUMBUS  home  page  is  at  http://www.itc.iinivie.ac.at/~hans/Coliimbus/coliimbus.html 
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Geometry  Optimizations 


Geometry  optimizations  were  performed  using  NWChem  version  4.0.1  on  both 
the  uranyl  and  uranate  ions  to  detennine  molecular  symmetry,  bond  lengths,  and 
vibrational  frequencies. 

All  geometry  optimizations  used  the  DRIVER  module  of  NWChem,  with  tight 
convergence  tolerance.  This  corresponded  to  a  root-mean-square  gradient  tolerance  of 
10‘5  and  a  maximum  gradient  of  1.5  •  1 0  5 .  Geometry  optimizations  were  performed 
without  symmetry  constraints,  in  order  to  determine  molecular  symmetry.  Subsequent 
geometry  optimizations  used  the  molecular  point-group  symmetry  determined  by  this 
geometry  optimization.  Only  those  geometries  which  yielded  all  positive  and  real 
vibrational  frequencies  at  both  DFT  and  MP2  levels  of  theory  were  considered  to  be  valid 
molecular  geometries. 

The  following  basis  sets  and  effective  core  potentials  were  used  for  the  uranium 
atom  during  the  geometry  optimizations  with  NWChem  4.0. 1 : 

•  Christiansen,  Ross,  and  Ennler  large-core  basis  and  ARECP  (Ennler  et  al, 
1991:  829-846) 

•  Kuchle,  Dolg,  Stoll  and  Preuss  small-core  basis  and  RECP  (Kuchle  et  al, 

1994:  7535-7542) 

For  the  oxygen  atoms,  the  following  basis  sets  and  core  potentials  were  used: 

•  Pacios  and  Christiansen  large-core  basis  and  RECP  (Pacios  et  al,  1985:  2664- 
2671) 
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•  Bergner,  Dolg,  Kiichle,  Stoll,  and  Preuss  small-core  basis  and  RECP  (Bergner 
etal,  1993:  1431-1441) 

•  Correlation-consistent,  polarized,  all-electron,  double-^  basis  set  (cc-pvdz) 
(Dunning,  1989:  1007-1023) 

•  Augmented,  correlation-consistent,  polarized,  double-^  (aug-cc-pvdz) 
(Kendall  etal,  1992:  6796-6806) 

Geometry  optimizations  were  performed  using  Hartree-Fock  SCF,  MP2,  and  DFT 
levels  of  theory.  For  DFT  optimizations  the  B3LYP  (Stephans,  1994:  1 1623),  Becke97 
(Becke,  1997),  and  Becke98  (Becke,  1998:  9624-9631)  exchange-correlation  functionals 
were  used.  Optimization  results  are  included  in  the  next  chapter. 

Unconstrained  uranyl  optimizations  resulted  in  a  linear  molecule.  As  a  result, 
subsequent  optimizations  were  constrained  using  D2h  and  D4h  point  group  symmetries. 

Unconstrained  uranate  optimizations  indicated  a  tetrahedral  and  distorted 
tetrahedral  geometries  for  the  gas  phase  ion. 

Overview  of  COLUMBUS  Calculations 

Once  an  appropriate  geometry  was  determined,  COLUMBUS  5.8.1  was  used  to 
compute  excited  electronic  state  energy  levels.  The  basic  procedure  for  COLUMBUS 
calculations  involved  running  a  series  of  FORTRAN  programs,  culminating  in  the  MR- 
CISD  calculation.  The  following  programs  were  run,  in  order,  to  compute  the  excited 
electronic  state  energy  levels:  ARGOS,  CNVRT,  SCFPQ,  MCDRT,  MCUFT,  MCSCF, 
MOFMT,  CIDRT,  CIUFT,  TRAN,  CISRT  and  CIUDG. 
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The  ARGOS  program  computes  the  one-  and  two-electron  atomic  orbital  integrals 
used  in  the  SCFPQ,  MCSCF,  and  CIUDG  programs.  The  CNVRT  program  converts 
these  integrals  into  PK  supermatrix  form,  to  speed  up  SCF  convergence.  The  SCFPQ 
program  performs  the  HF  self-consistent  field  calculations.  The  MCDRT  and  MCUFT 
programs  set  up  the  distinct  row  table  for  the  MCSCF  calculations.  The  MCSCF 
program  performs  the  multi-configuration  self-consistent  field  calculation,  while  the 
MOFMT  program  extracts  molecular  and  natural  orbital  coefficients  from  the  converged 
MCSCF  wave  function.  The  CIDRT  and  CIUFT  programs  set  up  the  distinct  row  table 
for  the  MR-CISD  calculation.  The  TRAN  program  converts  the  atomic  orbital  integrals 
into  molecular  orbital  integrals,  and  the  CISRT  program  sorts  these  integrals  into  groups 
of  zero-external,  one-extemal,  and  two-extemal  integrals  (Yabushita  et  al,  1999:  5797). 
The  CIUDG  program  then  diagonalizes  the  matrix-vector  products  of  the  Hamiltonian 
matrix  and  trial  vectors  in  order  to  obtain  the  ground  and  excited  electronic  state 
energies. 

COLUMBUS  Calculations  on  Uranyl  (U022+) 

Calculations  on  uranyl  were  performed  using  the  linear  D2h  molecular  point-group 
symmetry,  an  abelian  subgroup  of  Dooh-  Refer  to  Appendix  A  for  a  discussion  on  the 
symmetry  aspects  of  the  linear  uranyl  ion  and  a  definition  of  D^h  and  D2h  terms  and  states 
arising  from  electronic  configurations  in  uranyl. 

Because  of  symmetry,  the  uranyl  potential  energy  surface  (PES)  could  be 
characterized  by  a  single  parameter:  the  U-0  bond  length.  With  symmetry,  there  is  only 
one  symmetry  unique  oxygen  atom  in  the  uranyl  ion.  Varying  this  U-0  bond  length  was 
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equivalent  to  varying  the  symmetric  stretch  normal  mode  coordinate.  Calculations  were 
performed  using  both  the  large  and  small-core  uranium  RECPs.  The  small-core  RECP 
was  obtained  from  Dr.  Phil  Christiansen5.  The  small-core  basis  set  was  obtained  from 
Dr.  Russell  Pitzer’s  group  at  Ohio  State  University  (OSU).  The  small-core  basis  set  is  a 
generally  contracted  basis  set,  (7s8p7d5f4g)— >[7s6p7d2f4g]. 

Neither  the  small-core  RECP,  nor  its  associated  basis  set  have  been  published. 

The  only  difference  between  the  small-core  RECP  and  the  large-core  RECP  was  that  the 
small-core  potential  included  the  uranium  5d  shell  in  the  valence  space,  while  this  shell 
was  included  in  the  core  of  the  large-core  potential.  The  basis  set  used  for  both  the  large- 
and  small-core  uranyl  calculations  included  several  diffuse  p-functions  and  four 
polarization  g-functions. 

The  large  core  calculations  were  done  in  order  to  calibrate  the  calculations  and 
compare  results  to  those  obtained  by  OSU  (Zhang  et  al,  1999:  6880-6886).  The  small 
core  calculations  were  then  compared  to  the  large  core  uranyl  results,  and  both  results 
were  compared  to  uranyl  spectroscopic  data  (Denning,  1992:  217-275;  Rabinowitch  et 
al,  1964:  48).  This  comparison  allowed  for  qualitative  conclusions  to  be  drawn  on  the 
accuracies  of  the  large  core  RECP  versus  the  small  core  RECP. 

First,  a  series  of  SCF  calculations  were  performed  over  a  range  of  U-0  bond 
lengths.  The  minimum  to  this  potential  energy  surface  was  determined  from  a  4th  order 
polynomial  fit  to  these  SCF  energies.  Expanding  this  fitting  polynomial  in  a  Taylor 
series  about  the  interpolated  equilibrium  point,  then  computing  the  second  derivative 

5  Christiansen-style  RECPs  and  spin-orbit  potentials  can  be  downloaded  from 
http://www.clarkson.edu/~pac/reps.html 
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yielded  the  symmetric  stretch  vibrational  mode  force  constant,  from  which  the  symmetric 
stretch  vibrational  frequencies  were  computed.  This  fitting  procedure  was  used  to 
determine  the  potential  energy  surface  minima  and  symmetric  stretch  force  constants  in 
all  subsequent  calculations  on  uranyl  and  uranate. 

Following  this  series  of  SCF  ground  state  calculations,  a  series  of  MCSCF 
calculations  were  performed,  averaging  the  3 ex]  ground-state  with  the  3cr'  1 8]  and 

3cr'  1  (p'n  low  lying  excited  state  configurations  arising  from  singlet  and  triplet  Ag  and 
triplet  Big,  B2g,  and  B3g  symmetries.  The  2nu  and  3 nu  orbitals  were  thoroughly  mixed 
in  this  calculation,  in  a  manner  similar  to  that  done  in  the  original  calculation  at  OSU 
(Zhang  et  al,  1999:  6883).  The  uranium  5d  orbitals  were  frozen  in  the  small-core 
MCSCF  calculation.  Thus,  the  uranium  5d  orbitals  were  forced  to  be  occupied  in  all 
references,  however,  the  5d  electrons  were  included  in  the  correlation  calculations. 
Freezing  these  electrons  prevented  excitations  from  the  5d  shell.  Also  frozen  were  the 
highest  three  Ag  virtual  orbitals  and  the  highest  Big,  B2g,  and  B3g  virtual  orbitals.  This 
state-averaged  MCSCF  wave  function  over  the  ground-  and  first  low-lying  excited  states 
was  used  in  order  to  obtain  a  balanced  description  of  the  ground  and  low-lying  excited 
states.  Both  the  large-core  and  small-core  MCSCF  calculations  were  over  a  restricted 
active  space  consisting  of  two  electrons  in  five  orbitals,  where  the  first  orbital  was 
restricted  to  have  at  least  one  electron. 

Next,  the  converged  MCSCF  wave  function  was  used  as  the  starting  point  for  a 
MR-CISD  computation.  Spin-orbit  effects  were  included  in  the  calculation  during  this 
step.  The  same  three  references  used  in  the  MCSCF  calculation  were  also  used  in  the 
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MR-CISD  calculation,  and  the  same  core  and  virtual  orbitals  were  frozen.  The  MR- 
CISD  calculations  were  performed  using  both  molecular  and  natural  orbitals  extracted 
from  the  MCSCF  wave  function.  For  the  large-core  calculation,  the  MR-CISD  consisted 
of  roughly  3.5  million  CSFs.  The  small-core  calculation  had  nearly  9  million  CSFs. 

As  shown  in  the  Appendix  A,  all  states  arising  from  the  Eg,  ‘  Ag,  and  Og  states 
have  either  Ag  symmetry  (X„  states),  degenerate  Ag+Big  (Ag  states)  symmetry,  or 
degenerate  B2g+B3g  (Og)  symmetries.  In  the  co-co  spin-orbit  coupling  scheme6,  the  states 
with  Ag  correspond  to  (f,  states,  those  states  corresponding  with  degenerate  Ag+Big 

correspond  with  2  and  4  states,  while  degenerate  B2g+B3g  states  correspond  with  I ,, 
or  3g  states. 

Because  of  this  degeneracy,  the  lowest  six  excited  states  could  be  uniquely 
determined  by  examining  the  first  four  roots  with  Ag  symmetry  and  the  first  three  roots 
with  B2g  symmetry.  Degeneracy  of  these  roots  was  checked  by  computing  the  first  four 
Big  roots  at  the  SCF  minimum  energy  bond  length.  Using  the  minimum  MR-CISD 
energies  from  the  fourth-order  polynomial  fit  for  each  state,  the  electronic  transition 
energies  were  then  computed  and  compared  with  experimental  data.  State  assignments 
were  made  by  examining  the  Cl  wave  function  for  each  root.  Details  for  this  state 
assignment  procedure  are  contained  in  Appendix  C.  Results  from  these  calculations  are 
included  in  the  next  chapter. 


6  The  spin-orbit  coupling  scheme  where  electron  correlation  is  small  with  respect  to  spin-orbit  coupling  in 
atoms  is  j-j  coupling.  This  coupling  scheme  is  (0(0  -coupling  in  linear  molecules. 
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COLUMBUS  Calculations  on  Uranate  (U042  ) 

Geometry  optimizations  on  the  uranate  ion  using  both  MP2  and  DFT  indicated 
were  dependent  on  the  geometry  of  the  molecule.  Optimizations  beginning  from  a  square 
planar  molecule  converged  to  a  square  planar  molecule.  Optimizations  that  began  with  a 
tetrahedral  geometry  converged  to  a  tetrahedral  geometry.  Optimizing  the  molecular 
geometry  using  an  asymmetric  starting  point  resulted  in  a  distorted  tetrahedral  molecule 
using  both  DFT  and  MP2.  The  square  planar  and  tetrahedral  geometries  are  likely 
stationary  points  on  the  potential  energy  surface.  The  tetrahedral  geometry  was  used  in 
subsequent  calculations,  because  the  potential  energy  surface  could  be  characterized  by  a 
single  parameter. 

Either  D2  or  C-2V  point- groups  would  be  appropriate  abelian  subgroups  of  the  Td 
symmetry  point-group,  and  the  calculation  on  uranate  was  carried  out  using  D2  point- 
group  symmetry.  See  Appendix  B  for  a  discussion  symmetry  aspects  of  the  tetrahedral 
uranate  ion  and  a  definition  of  terms  and  states  arising  from  the  Td  and  D2  point-groups. 

As  was  the  case  with  the  uranyl  ion,  the  uranate  ion  potential  energy  surface 
(PES)  could  be  characterized  by  a  single  U-0  bond  length  in  Td  and  D2  point-group 
symmetry.  Since  there  is  only  one  symmetry  unique  oxygen  atom  in  a  tetrahedral 
geometry,  varying  this  bond  length  simultaneously  varied  the  bond  lengths  of  the  other 
oxygen  atoms.  This  procedure  is  equivalent  to  varying  the  symmetric  stretch  nonnal 
mode  coordinate  in  the  molecule.  Varying  one  U-0  bond  length  varied  all  U-0  bonds 
simultaneously  by  equal  amounts.  Based  on  the  spectroscopic  accuracy  obtained  from 
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the  small-core  uranyl  calculation,  the  uranate  calculations  were  performed  using  only  the 
small-core  uranium  RECP. 

The  initial  ARGOS  input  was  generated  automatically  using  the  IARGOS 
program.  This  method  resulted  in  several  errors.  First,  the  IARGOS  program  generated 
symmetry  adapted  linear  combinations  of  the  f-  and  g-  functions  appropriate  for  axially 
symmetric  molecules.  For  the  tetrahedral  uranate  molecule,  a  symmetry  adapted  linear 
combination  using  cubic  f-  and  g-  functions  was  appropriate.  Additionally,  the  program 
created  an  incorrect  symmetry  operation  generator  for  the  molecule.  These  errors  were 
subtle,  and  were  difficult  to  detect.  They  also  resulted  in  electronic  energies  that  did  not 
exhibit  the  appropriate  symmetry  appropriate  for  a  tetrahedral  molecule  in  D2. 

Another  problem  was  identified  with  the  basis  set.  The  same  basis  set  used  for 
the  uranyl  calculations  was  used  for  the  uranate  calculations.  SCF  and  MCSCF 
calculations  resulted  in  exaggerated  uranium  7p  population,  which  was  unexpected. 
Based  on  the  uranyl  results,  the  uranium  5f  orbitals  were  expected  to  play  a  large  role  in 
the  excited  states  of  uranate.  Subsequent  analysis  revealed  that  the  large  number  of 
diffuse  p-functions  and  g-polarization  functions  biased  the  calculation  in  favor  of  the 
uranium  7p  orbitals.  The  problem  was  fixed  by  removing  the  diffuse  p-functions  and 
removing  all  but  one  g-polarization  function.  The  absence  of  non-bonding  orbitals  in  the 
tetrahedral  uranate  revealed  the  impact  of  the  diffuse  p-  and  polarization  g-  functions. 
Similar  effects  did  not  show  up  in  the  uranyl  calculation,  because  of  the  presence  of  the 
non-bonding  15u  and  l(pu  orbitals.  In  the  uranyl  calculation,  the  basis  set  impacted  those 
anti-bonding  states  that  were  much  higher  in  energy  than  the  low  lying  states  of  interest. 
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Again,  a  series  of  SCF  calculations  were  performed  over  a  range  of  symmetric 
stretch  nonnal  mode  coordinate  U-0  bond  lengths.  Fourth-order  polynomial  fits  were 
used  to  find  the  minimum  energies  and  bond  lengths  for  SCF,  MCSCF,  and  MR-CISD 
data.  Expanding  this  fitting  polynomial  in  a  Taylor  series  about  the  equilibrium  point, 
then  finding  the  second  derivative  yielded  the  symmetric  stretch  force  constants  for  this 
nonnal  mode,  which  in  turn  yielded  the  symmetric  stretch  vibration  frequencies.  From 
the  SCF  calculation,  the  uranate  HOMO  proved  to  be  the  5t2  orbital,  while  the  LUMO 
was  the  6t2  orbital.  The  SCF  results  showed  a  significant  HOMO-LUMO  energy  gap, 
unlike  uranyl.  This  was  due  to  the  absence  of  non-bonding  orbitals,  unlike  the  uranyl 
molecule.  The  large  HOMO/LUMO  energy  gap  in  the  uranate  molecule  indicated  a 
single-reference  calculation  would  suffice,  as  opposed  to  uranyl,  where  a  state-averaged, 
multi-reference  description  of  the  ground  state  was  more  appropriate. 

Following  this  series  of  SCF  ground  state  calculations,  a  series  of  MCSCF 
calculations  were  performed,  using  the  ground  state  reference  only.  In  this  aspect,  the 
MCSCF  calculation  was  identical  to  the  HF-SCF  calculation.  The  only  reason  for 
performing  this  calculation  was  to  obtain  natural  orbitals  for  use  in  the  MR-CISD 
calculation.  As  was  done  in  the  uranyl  calculations,  the  uranium  5d  orbitals  were  frozen 
in  this  calculation,  meaning  they  were  occupied  in  all  references.  While  excitations  from 
the  uranium  5d  shell  were  not  allowed,  the  electrons  were  allowed  to  participate  in  the 
correlation  calculation. 

Using  the  natural  orbital  coefficients  from  the  MCSCF  calculation  as  a  starting 
point,  a  single -point  MR-CISD  calculation  was  performed  on  the  uranate  ion  at  the 
interpolated  SCF  minimum  U-0  bond  length.  The  references  used  in  the  uranate 
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calculation  were  similar  in  nature  to  the  uranyl  calculation.  Three  references  were 
chosen:  the  ground  state,  and  those  references  consisting  of  single  excitations  into  the 
low-lying,  unoccupied  uranium  5f  orbitals  (8ai,  6t2),  with  the  requirement  of  five 
electrons  remaining  in  the  HOMO.  A  small  MR-CISD  calculation  with  all  occupied 
orbitals  except  for  the  5t2  HOMO  frozen,  and  48  virtual  orbitals  frozen  was  performed  in 
order  to  characterize  the  low-lying  excited  states.  This  small  calculation  yielded  the 
symmetry  expected  of  a  tetrahedral  molecule,  and  gave  a  qualitative  picture  of  the  nature 
and  ordering  of  the  low-lying  excited  states. 

It  was  shown  in  Appendix  B  that  a  large  number  of  possible  states  arise  from  the 
electron  configurations  5 1\  ,  5 1\  6 1\ ,  and  5 ts2l  6 1\ .  Unlike  the  calculation  with  uranyl,  the 

nature  and  ordering  of  the  low  lying  excited  states  for  uranate  were  unknown.  Appendix 
C  contains  the  details  for  the  double-group  state  assignment  for  the  uranate  ion  ground 
and  excited  electronic  states. 
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IV.  Results  and  Discussion 


In  this  chapter,  the  results  of  the  geometry  optimizations  using  NWChem  4.0. 1,  as 
well  as  the  results  of  the  COLUMBUS  5.8.1  SCF,  MCSCF,  and  MR-CISD  calculations 
on  both  the  uranyl  and  uranate  ions  will  be  presented. 

Uranyl  Geometry  Optimization  Results 

The  initial  results  of  the  geometry  optimizations  for  the  uranyl  ion  are  presented  in  Table 
1 .  Also  listed  is  a  high-fidelity  theoretical,  four-component,  all-electron  calculation 
incorporating  electron  correlation  at  roughly  a  triple  and  quadruple  excitation  level.  This 
calculation  represents  the  current  state-of-the-art  in  relativistic  quantum  chemistry 
calculations  on  uranyl  (de  Jong  et  al,  1999:  41-52);  unfortunately,  excited  states  were  not 
investigated  in  the  study  by  de  Jong  et  al.  Table  2  and  Table  3  lists  the  interpolated 
minimum  bond-length  and  symmetric  stretch  vibrational  frequency  computed  from  the 
COLUMBUS  large-  and  small-core  potential  energy  surfaces  for  comparison  with  the 
NWChem  4.0.1  results.  Several  experimental  results  measuring  the  vibrational 
frequencies  of  uranyl  are  also  listed  for  comparison  (Toth  et  al,  1981:  547-549;  Denning 
et  al,  1992:  216-275).  The  first  experimental  results  from  Toth  et  al  are  in  an  aqueous 
HNO3  environment,  while  the  second  experiment  was  perfonned  on  a  crystal  of 
CS2UO2CI4.  It  is  known  that  the  uranyl  vibrational  frequencies  can  be  perturbed  by  their 
crystalline  or  aqueous  environment.  However,  the  strength  of  the  uranyl  bond  minimizes 
the  impact  of  crystal  field  perturbations  on  the  uranyl  electronic  spectra.  Symmetric 
stretch  vibrational  frequencies  vary  from  887.9  cm"1  in  RbUCLCNCfifi  to  808  cm'1  in 
KMJCLCCCL^,  while  the  uranyl  fluorescent  series  ranges  from  21 199.4  cm'1  in 
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RbU02(N03)3  to  20096.3  cm'1  in  CS2UO2CI4  (Rabinowitch  et  al,  1964:  48).  The  largest 
crystal  field  perturbation  on  the  symmetric  stretch  vibrational  frequency  for  uranyl  is 
roughly  80  cm'1,  while  the  maximum  crystal  field  splitting  of  the  electronic  fluorescent 
spectra  appears  to  be  roughly  1 103  cm'1. 
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Table  1.  Uranyl  Geometry  Optimization  Results 


Uranium 

RECP/basis 

Oxygen  basis 

method 

Constraint 

Symmetry 

bond  length  (A) 

symmetric 
stretch  tcm'1) 

Kiichle  et  al, 
1994 

aug-cc-pvdz 
(Kendall  et  al, 
1992) 

DFT 

B3LYP 

D’h 

1.6997 

1037 

Kiichle  et  al, 
1994 

aug-cc-pvdz 
(Kendall  et  al, 
1992) 

DFT 

B3LYP 

fKii 

1.6997 

1037 

Kiichle  et  al, 
1994 

aug-cc-pvdz 
(Kendall  et  al, 
1992) 

DFT  B3LYP 

Ci 

1.6487 

1219 

Kiichle  et  al, 
1994 

aug-cc-pvdz 
(Kendall  et  al, 
1992) 

DFT  Becke97 

D’h 

1.6950 

1051 

Kiichle  et  al, 
1994 

aug-cc-pvdz 
(Kendall  et  al, 
1992) 

DFT  Becke98 

D’h 

1.6924 

1059 

Kiichle  et  al, 
1994 

Bergner  et  al, 
1993 

MP2 

Cl 

1.7758 

906 

Ennler  et  al, 
1991 

Pacios  et  al,  1985 

DFT  B3LYP 

Ci 

1.7257 

921 

Ennler  et  al, 

1991 

Pacios  et  al,  1985 

DFT  B3LYP 

Ci 

1.7040 

991 

Ennler  et  al, 
1991 

aug-cc-pvdz 
(Kendall  et  al, 
1992) 

DFT  B3LYP 

D411 

1.6850 

1026 

Ennler  et  al, 

1991 

aug-cc-pvdz 
(Kendall  et  al, 
1992) 

DFT  Becke97 

D’h 

1.6834 

1038 

Ennler  et  al, 
1991 

aug-cc-pvdz 
(Kendall  et  al, 
1992) 

DFT  Becke98 

D’h 

1.6810 

1045 

Ennler  et  al, 
1991 

Pacios  et  al,  1985 

MP2 

Cl 

1.7572 

891 

Ennler  et  al, 
1991 

aug-cc-pvdz 
(Kendall  et  al, 
1992) 

HF  SCF 

D’h 

1.6356 

1097 

Ennler  et  al, 

1991 

aug-cc-pvdz 
(Kendall  et  al, 
1992) 

HF  SCF 

Cl 

1.6356 

1099 

Ennler  et  al, 
1991 

aug-cc-pvdz 
(Kendall  et  al, 
1992) 

HF  SCF 

D’h 

1.6341 

1097 

4-component  all¬ 
electron 
theoretical 
calculation7 

DHF  +  CCSD(T) 

D2h2  double 
group 

1.715 

974 

7  de  Jong  et  al,  1999 
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Table  2.  COLUMBUS  Large-Core  Uranyl  Ground-State  Calculation  Results 


Method 

bond  length 

(A) 

symmetric 
stretch  (cm'1) 

HF  SCF 

1.6465 

996 

MCSCF 

1.6522 

1065 

MR-C1SD 

1.6679 

1038 

Table  3.  COLUMBUS  Small-Core  Uranyl  Ground-State  Calculation  Results 


Method 

bond  length 

(A) 

symmetric 
stretch  (cm'1) 

HF  SCF 

1.6625 

999 

MCSCF 

1.6717 

1087 

MR-C1SD 

1.6869 

1062 

2 

Table  4  shows  the  geometry  optimization  results  for  the  UO4  '  ion.  Both  DFT  and 
MP2  geometry  optimizations  yielded  tetrahedral  (Td)  when  started  in  a  tetrahedral 
geometry.  Both  MP2  and  DFT  yielded  square  planar  molecules  when  started  from  a 
square  planar  geometry.  An  interesting  result  occurred  when  MP2  and  DFT 
optimizations  started  with  an  asymmetric  geometry.  The  DFT  optimizations  converged 
to  a  very  flat  Did  geometry,  almost  square  planar.  The  potential  energy  surface  of  this 
geometry  was  very  flat,  and  the  geometry  optimization  could  not  converge  to  the 
specified  gradient  tolerance  of  1.0  •  10-'  . 
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Table  4.  Uranate  NWChem  4.0. 1  Geometry  Optimization  Results 


U  basis 

O  basis 

method 

Svmmetrv 

bond  length 

(A) 

symmetric 
stretch  (cm'1! 

Ermler  et  al, 
1991 

aug-cc-pvdz 
(Kendall  et  al, 
1992) 

MP2 

Ci 

2.0164 

tetrahedral 

N/A 

Ermler  et  al, 
1991 

cc-pvdz 
(Dunning  et 
al,  1989) 

DFT  B3LYP 

Ci 

1.9764 

tetrahedral 

724 

Ermler  et  al, 
1991 

Pacios  et  al, 
1985 

DFT  B3LYP 

Cl 

1.9918 

tetrahedral 

707 

Ermler  et  al, 
1991 

aug-cc-pvdz 
(Kendall  et  al, 
1992) 

DFT  B3LYP 

Cl 

1.9851 

tetrahedral 

713 

Kiichle  et  al, 
1994 

aug-cc-pvdz 
(Kendall  et  al, 
1992) 

DFT  B3LYP 

Td 

1.9676 

750 

Ermler  et  al, 
1991 

aug-cc-pvdz 
(Kendall  et  al, 
1992) 

DFT  B3LYP 

Cl 

1.933 

nearly  square 
planar 

N/A 

Large-core  Uranyl  (U022+)  Results 

COLUMBUS  5.8.1  large-core  calculation  results  are  presented  next.  Figure  4 
shows  the  results  of  the  FIF  SCF  calculations  as  a  function  of  bond  length.  Because  of 
symmetry,  this  bond  length  represents  the  symmetric  stretch  normal  mode  coordinate  of 
the  molecule.  Varying  this  bond  length  simultaneously  varied  both  oxygen-uranium 
bond  lengths  by  equal  amounts. 
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Large  Core  SCF  PES 


Bond  Length  (A) 


Figure  4.  Large-core  Uranyl  SCF  Potential  Energy  Surface 
A  fourth-order  polynomial  fit  to  this  data  yielded  a  minimum  total  electronic 
energy  at  a  bond  length  of  1.647  A.  The  second  derivative  of  the  Taylor  series  expansion 
about  the  equilibrium  point  yielded  the  force  constant  for  the  symmetric  stretch  normal 
mode.  Table  5  contains  the  symmetric  stretch  vibrational  frequency,  calculated  from  this 
force  constant.  The  frequency  was  scaled  by  90%,  as  is  frequently  appropriate  for  SCF 
frequency  calculations  (Levine,  2000:  703-704). 

Table  5.  Large-core  Uranyl  SCL  Symmetric  Stretch  Vibrational  Frequencies 


State 

Equilibrium  Bond 

Length  (A) 

Scaled  Vibrational 

Frequency  (cm1) 

1.6465 

996 

64 


Next,  the  state-averaged  MCSCF  calculation  over  the  3crJ  ,  3cr'  Iff ,  and  3cr'  \(p[ 


electron  configuration  references  are  shown  in  Figure  5. 


Large  Core  MCSCF  PES 
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Figure  5.  Large-core  Uranyl  MCSCF  Potential  Energy  Surface 
Again,  equilibrium  bond  lengths  were  obtained  from  fourth-order  polynomial  fits 
to  the  MCSCF  total  electronic  energy  data,  while  the  symmetric  stretch  vibrational 
frequencies  were  obtained  from  the  second  derivative  of  a  Taylor  series  expansion  about 
the  equilibrium  point.  Table  6  lists  the  equilibrium  bond  lengths  and  symmetric  stretch 
vibrational  frequencies  obtained  from  the  MCSCF  data. 
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Table  6.  Large-core  Uranyl  MCSCF  Symmetric  Stretch  Vibrational  Frequencies 


State 

Equilibrium  Bond 
Length  (A) 

Symmetric  Stretch 
Vibrational 
Frequency  (cm'1) 

1.6521 

1065 

3\ 

1.7082 

861 

3®, 

1.7310 

815 

The  MR-CISD  calculation  was  performed  using  the  same  three  references  used  in 
the  MCSCF  calculation.  Both  molecular  and  natural  orbital  coefficients  were  extracted 
from  the  MCSCF  wave  function  and  used  as  the  initial  guess  for  the  CISD  calculation. 

Figure  6  shows  the  MR-CISD  potential  energy  surface  obtained  using  the 
molecular  orbital  coefficients  as  the  initial  guess  to  the  MR-CISD  wave  function. 
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Large-core  Uranyl  MR-CISD  Electronic  States 
from  MCSCF  Molecular  Orbitals 
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Figure  6.  Large-core  MR-CISD  Uranyl  Electronic  States  from  Molecular  Orbitals 

Figure  7  shows  the  MR-CISD  potential  energy  surface  obtained  by  using  the 
natural  orbital  coefficients  as  the  initial  guess  to  the  MR-CISD  wave  function.  As 
explained  in  the  theory  section,  Cl  expansions  using  natural  orbitals  are  generally  more 
accurate  then  those  obtained  from  molecular  orbitals.  The  improved  accuracy  of  the 
natural  orbital  Cl  expansion  is  illustrated  by  comparing  the  shapes  of  the  MR-CISD 
potential  energy  surfaces  obtained  using  molecular  and  natural  orbital  starting  guesses. 
The  potential  energy  surfaces  obtained  using  natural  orbital  expansions  yielded  smoother 
potential  energy  surfaces,  while  the  molecular  orbital  Cl  expansion  resulted  in  distorted 
potential  energy  surfaces  for  the  excited  states. 
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Large-core  MR-CISD  Electronic  States  from 
MCSCF  Natural  Orbitals 
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Figure  7.  Large-core  MR-CISD  Uranyl  Electronic  States  from  Natural  Orbitals 

Equilibrium  bond  lengths  and  symmetric  stretch  vibrational  frequencies  for  the 
large-core  MR-CISD  electronic  states  obtained  from  natural  orbitals  are  listed  in  Table  7. 
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Table  7.  Large-core  Uranyl  MR-CISD  Symmetric  Stretch  Vibrational  Frequencies 


State 

Equilibrium  Bond 
Length  (A) 

Symmetric  Stretch 
Vibrational 
Frequency  (cm'1) 

1.6630 

1038 

1.6967 

840 

1.7022 

804 

3, 

1.7065 

802 

2g 

1.7288 

839 

3g 

1.7252 

839 

4g 

1.7483 

809 

Finally,  the  composition  of  the  MR-CISD  wave  function  fonned  from  natural 
orbitals  is  listed  in  Table  8.  The  effect  of  the  spin-orbit  potential  in  the  Cl  calculation  is 
apparent  here.  Compare  the  number  of  states  arising  from  the  MCSCF  and  MR-CISD 
calculations  spin-orbit.  The  MR-SOCI  calculation  splits  the  MCSCF  states  in  Figure  5, 
resulting  in  those  states  shown  in  Figure  7.  Compositions  do  not  necessary  sum  to  100%, 
due  to  a  large  number  of  contributions  from  configuration  state  functions  with  small 
coupling  coefficients.  Additionally,  coupling  coefficients  alternated  in  sign.  Only  those 
states  with  expansion  coefficients  larger  than  0.1  are  listed  in  Table  8. 
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Table  8.  Large-core  Uranyl  Wave  Function  Compositions 


State 

Large-Core  Wave 
Function  Composition 

Zhang  et  al 
(Zhang,  1999:  6884) 

O 

oq  + 

82%  'z; 

83%  ‘z; 

84%  3Ag 

84%  3Ag 

58%  3Ag 

62%  3Ag 

26%  3©g 

22%  3©g 

3, 

52%  3Ag 

55%  3Ag 

26%  3©g 

23%  3©g 

6%  3©g 

6%  1©g 

2, 

25%  3Ag 

21%  3Ag 

3%  ‘Ag 

3%  'Ag 

56%  3©g 

60%  3©g 

3, 

32%  3Ag 

28%  3Ag 

48%  3©g 

51%  3©g 

4%  3©  g 

4%  1©g 

4 

84%  3<f> 

84%  3®, 

g 

g 

g 

Finally,  the  adiabatic  electronic  transition  energies  were  computed  from  the  large- 
core  MR-CISD  equilibrium  energies.  These  results  were  compared  to  experimental 
results  obtained  by  Denning  et  al  on  CS2UO2CI4  (Denning,  1992)  and  CsU02(N03)3 
(Denning  et  al,  1979).  Shown  in  Table  9,  it  is  apparent  that  the  MR-CISD  method  using 
the  large-core  RECP  models  the  optical  electronic  spectra  of  uranyl  rather  well,  to  within 
about  3%  of  the  experimental  data,  on  average. 
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Table  9.  Electronic  Transition  Energies  from  Large-core  MR-CISD  Results 


Electronic 

state 

Equilibrium 
bond  length  (A) 

adiabatic 

AE  (cm1) 

CsjUOzCk 

adiabatic  % 
difference 

og+ 

1.6679 

ls 

1.7324 

21430 

20096 

3.2% 

2g 

1.7376 

22151 

20861 

3.0% 

3g 

1.7402 

23378 

22051 

2.9% 

2g 

1.7471 

24637 

22578 

4.4% 

3S 

1.7458 

26872 

26222 

1.2% 

4g 

1.7547 

28736 

27738 

1.8% 

average 

2.9% 

Table  10.  Large-core  Elranyl  Results  from  Zhang  et  al  (Zhang,  1999:  6884) 


Electronic 

state 

Equilibrium  bond  length 

(A) 

Symmetric  Stretch 
Vibrational  Frequency 
(cm1) 

adiabatic  AE 
(cm1) 

0g+ 

1.668 

1103 

lg 

1.733 

867 

20719 

2g 

1.739 

845 

21421 

3g 

1.742 

847 

22628 

2g 

1.749 

900 

23902 

3g 

1.747 

898 

26118 

_ 4g _ 

1.755 

880 

27983 
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The  comparison  between  this  calculation  and  the  one  performed  by  Zhang  et  al 
(Zhang,  1999:  6884)  is  quite  favorable.  Bond  lengths  match  closely,  as  do  symmetric 
stretch  vibrational  frequencies.  Adiabatic  transition  energies  are  slightly  different; 
however,  this  difference  is  not  unexpected  because  of  the  slightly  different  treatment  of 
the  ground  and  excited  states  versus  the  calculation  performed  by  Zhang  et  al. 


Small-core  Uranyl  (UC>22+)  Results 


The  small-core  HF  SCF  potential  energy  surface  is  displayed  in  Figure  8. 


Small  Core  SCF  PES 
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Figure  8.  Small-core  Uranyl  SCF  Potential  Energy  Surface 
As  was  done  in  the  large-core  calculations,  equilibrium  bond  lengths  and 
minimum  energies  were  computed  from  fourth-order  data  fits,  while  symmetric  stretch 
vibrational  frequencies  were  computed  from  second  derivative  of  the  Taylor  series 
expansion  about  the  equilibrium  point.  Again,  the  symmetric  stretch  vibrational 
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frequency  from  the  HF  SCF  potential  energy  surface  was  scaled  by  90%.  This  data  is 


listed  in  Table  11. 


Table  11.  Small-core  Uranyl  SCF  Equilibrium  Bond  Length  and  Symmetric  Stretch 

Vibrational  Frequency 


State 

Equilibrium  Bond 

Length  (A) 

Scaled  Symmetric 

Stretch  Vibrational 

Frequency  (cm'1) 

1  y+ 

g 

1.6625 

999 

The  potential  energy  surfaces  obtained  from  the  small-core  MCSCF  calculation 
are  shown  in  Figure  9. 
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Figure  9.  Small-core  Uranyl  MCSCF  Potential  Energy  Surface 
Table  12  lists  the  equilibrium  bond  lengths  and  symmetric  stretch  vibrational 
frequencies  obtained  from  the  state-averaged  MCSCF  results. 
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Table  12.  Small-core  Uranyl  MCSCF  Symmetric  Stretch  Vibrational  Frequencies 


State 

Equilibrium  Bond 

Length  (A) 

Symmetric  Stretch 

Vibrational 
Frequency  (cm'1) 

1.6717 

1087 

1.7263 

903 

1.7442 

871 

As  was  done  with  the  large-core  calculation,  the  small-core  MR-CISD  calculation 
was  perfonned  using  both  molecular  and  natural  orbitals  as  starting  points.  The  results  of 
the  Cl  expansion  from  the  molecular  orbital  coefficients  is  shown  in  Figure  10.  In  this 
calculation,  the  molecular  orbital  coefficients  are  particularly  poor,  as  they  fail  to  yield 
bound  excited  states  for  uranyl.  This  situation  does  not  represent  reality. 


Small-core  Uranyl  MR-CISD  Electronic  States 
from  MCSCF  Molecular  Orbitals 
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Figure  10.  Small-core  Uranyl  MR-CISD  Electronic  States  From  Molecular  Orbitals 
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The  results  Cl  expansion  using  the  natural  orbitals  are  shown  in  Figure  11.  The 
difference  in  the  accuracies  between  molecular  and  natural  orbital  initial  guesses  had  a 
substantial  impact  on  the  final  MR-CISD  wave  function.  With  the  small-core 
calculation,  the  accuracy  difference  makes  the  difference  between  bound  and  unbound 
states. 

Small-core  Uranyl  MR-CISD  Electronic  States 
from  MCSCF  Natural  Orbitals 
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Figure  11.  Small-core  Uranyl  MR-CISD  Electronic  States  From  Natural  Orbitals 

The  equilibrium  bond  lengths  and  symmetric  stretch  frequencies  for  the  MR- 
CISD  electronic  states  are  listed  in  Table  13. 
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Table  13.  Small-core  Uranyl  MR-CISD  Symmetric  Stretch  Vibrational  Frequencies 


State 

Equilibrium  Bond 
Length  (A) 

Symmetric  Stretch 
Vibrational  Frequency  (cm'1) 

o; 

1.6869 

1062 

1.7470 

882 

2S 

1.7571 

845 

3, 

1.7592 

844 

2g 

1.7553 

885 

3g 

1.7546 

883 

4g 

1.7649 

856 

Next,  the  small-core  uranyl  MR-CISD  wave  function  composition  is  displayed  in 
Table  14.  The  large-core  uranyl  wave  function  composition  for  the  same  state  is  also 
listed  for  comparison.  The  only  fundamental  difference  between  the  wave  functions 
occurs  in  the  2g  states. 
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Table  14.  Small-  and  Large-core  Uranyl  Wave  Function  Compositions 


State 

Small-Core  MR-CISD 
Wave  Function 
Composition 

Laree-Core  MR-CISD 
Wave  Function 
Composition 

83%  lZ+g 

82%  'A; 

84%  3Ag 

84%  3A? 

2g 

37%  3Ag 

58%  3A? 

46%  3©g 

26%  3©g 

3g 

40%  3Ag 

52%  3A? 

38%  3©g 

26%  3©g 

10%  'cD? 

6%  lOg 

2g 

44%  3Ag 

25%  3Ag 

4%  'A? 

3%  'Ag 

35%  3©g 

56%  3©g 

47%  3Ag 

32%  3Ag 

35%  3©g 

48%  3©g 

2%  lOg 

4%  !©g 

4 

84%  3© 

84%  3© 

g 

g 

g 

Finally,  Table  15  shows  the  adiabatic  excited  state  transition  energies  computed 
from  the  MR-CISD  equilibrium  energies  and  compared  with  the  CS2UO2CI4  (Denning, 
1992)  and  CsU02(N03)3  (Denning  et  al,  1979)  experimental  data. 
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Table  15.  Electronic  Transition  Energies  from  Small-core  MR-CISD  Results 


State 

Equilibrium 

Bond 

Length  (A) 

adiabatic 
AE  (cm'1) 

Cs7.UO7.Cl4 

adiabatic 
%  diff 

°g+ 

1.6869 

u 

1.7470 

20813 

20096 

1.8% 

2g 

1.7571 

20951 

20861 

0.2% 

38 

1.7592 

22189 

22051 

0.3% 

2g 

1.7553 

23549 

22578 

2.1% 

3g 

1.7546 

26251 

26222 

0.1% 

4g 

1.7649 

27729 

27738 

0.0% 

average  AE 

0.7% 
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Table  16.  Comparison  Between  Theoretical  and  Experimental  Uranyl  Fluorescent 
Electronic  Spectra  and  Symmetric  Stretch  Vibrational  Frequencies  (Rabinowitch  et  al, 

1964:  48) 


Compound 

Fluorescent  Series 

(cm1) 

Symmetric  stretch  vibrational 

frequency  (cm1) 

Rb2U02Cl4  2H20 

19961 

831.8 

K2UO2CI4  2H20 

19970 

831.6 

CS2UO2CI4 

20096 

836.1 

K2U02(S04)2 

20389 

827 

Rb2U02(S04)2 

20390 

838.2 

U02(N03)2  6H20 

20578 

863.9 

Cs2U02(S04)2  3H20 

20594 

860.5 

Large  core  (OSU) 

20719 

867 

U02(N03)2  3H20 

20779 

874.0 

Rb2U02(N03)4 

20808 

887.9 

Small  core  (this  work) 

20813 

882 

K2U02(N03)4 

20818 

870.3 

K2U02(C03)2 

20943 

808.0 

PbU02(CH3C00)4 

20958 

853.0 

CsU02(CH3C00)3 

20992 

842.2 

RbU02(CH3C00)3 

21049 

852.1 

NH4U02(CH3C00)3 

21056 

847.0 

CsU02(N03)3 

21090 

884.0 

NH4U02(N03)4 

21098 

885.7 

NaU02(CH3C00)3 

21135 

855.2 

KU02(N03)3 

21183 

875.5 

RbU02(N03)3 

21199 

887.9 

Large  core  (this  work) 

21430 

840 

(NH4)2U02(S04)2 

23358 

840 

Experiment  average 

20884 

854.9 

The  small-core  uranyl  MR-CISD  energies  agrees  more  closely  with  the 
CS2UO2CI4  experimental  data  than  the  large-core  results,  coming  within  less  than  1%  of 
the  experimental  values  on  average.  Additionally,  the  experimental  transition  energies 
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are  component  averaged,  due  to  the  crystal  field  splitting  observed  in  the  experimental 
spectroscopy. 

Based  on  the  orbital  coefficients  obtained  from  the  SCF  wave  function,  the  3 cr2 

LUMO  consists  of  primarily  uranium  5d,  and  7s  orbitals  bonding  with  oxygen  Is  orbitals. 
There  is  also  a  small  oxygen  2p  component.  The  first  excited  states,  arising  from  the 
3cr'  \S'U  electron  configuration,  consist  almost  entirely  of  uranium  5f  orbitals,  indicating  a 
charge  transfer  process  between  the  axial  oxygen  atoms  and  the  uranium  atom.  The  last 
three  excited  states,  arising  from  the  3cr'  1  cp\  electron  configuration  also  are  substantially 

of  uranium  5f  character.  Thus,  the  lowest  lying  excited  states  of  uranyl  all  involve 
charge  transfer  from  oxygen  is  and  2p  orbitals  to  uranium  5f  orbitals. 

Small-core  Uranate  (U042')  Results 

Based  on  the  accuracy  of  the  small-core  uranyl  results  comparison  with  the 
precise  CS2UO2CI4  experimental  data,  all  uranate  ion  calculations  were  performed  using 
only  the  small-core  RECP.  The  SCF  equilibrium  bond  length  was  computed  from  fourth- 
order  polynomial  fits  to  the  raw  data,  while  the  symmetric  stretch  (breathing-mode) 
vibrational  frequencies  were  obtained  from  the  second  derivative  of  the  Taylor  series 
expansion  about  the  equilibrium  point.  The  resulting  HF  SCF  minimum  bond  length  was 
used  for  the  single-point  MR-CISD  energy  calculation.  Figure  12  shows  the  HF  SCF 
potential  energy  surface  computed. 
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Small-core  U042'  SCF  Potential  Energy  Surface 
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Figure  12.  Uranate  SCF  Potential  Energy  Surface 

The  uranate  ground  state,  since  it  is  a  closed  shell  molecule,  is  the  completely 
symmetric  singlet  state,  or  1  Ai  in  Td  symmetry.  The  SCF  breathing  frequency  was  scaled 
by  90%.  SCF  results  for  the  equilibrium  bond  length  and  breathing  frequency  are  shown 
in  Table  17. 
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Table  17.  Small-core  Uranate  SCF  Symmetric  Stretch  Vibrational  Frequencies 


State 

Equilibrium  Bond 
Length  (A) 

Symmetric  Stretch 
Scaled  Vibrational 
Frequency  (cm'1) 

lA 

1.9575 

780 

The  small-core  uranate  MCSCF  computation  used  only  the  ground  state  as  a 
reference.  The  decision  to  use  only  a  single  reference  was  based  upon  the  relatively  large 
HOMO-LUMO  energy  gap  found  in  the  SCF  calculation. 

Because  the  MCSCF  calculation  only  used  the  ground  state  as  a  reference,  its 
results  should  mirror  the  HF  SCF  calculation.  In  contrast,  the  uranyl  MCSCF  calculation 
was  state-averaged  over  the  ground  state  and  first  excited  states.  This  was  due  to  the 
relatively  small  separation  in  energy  between  the  uranyl  HOMO  and  LUMO.  The  only 
benefits  from  the  uranate  MCSCF  calculation  were  improved  virtual  orbital  energies  as 
opposed  to  the  HF  SCF  virtual  orbital  energies,  and  the  ability  to  extract  natural  orbitals 
from  the  MCSCF  wave  function  for  use  as  a  starting  wave  function  in  the  MR-CISD 
calculation. 

Using  the  natural  orbitals  extracted  from  the  single-reference  MCSCF  wave 
function  as  a  starting  point,  a  small  MR-CISD  calculation  was  performed  using  over 
several  references.  The  references  used  in  the  MR-CISD  calculation  included  single 
excitations  from  the  5 1\  HOMO,  which  was  mainly  an  oxygen  2p  orbital,  into  the  low- 
lying  uranium  5f  orbitals,  a  situation  analogous  to  the  uranyl  calculation.  These  low- 
lying  orbitals  were  the  8a}  and  6t2  orbitals.  The  lowest  48  virtual  orbitals  were  frozen  in 
this  small  Cl  calculation.  In  Appendix  B,  it  was  detennined  that  there  are  a  possible  13 
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states  of  A  symmetry,  including  the  ground  state,  and  a  total  of  12  possible  states  of  Bi, 
Bo,  and  B3  symmetry  arising  from  the  three  references  used.  The  number  of  CSFs  of  A 
symmetry  in  D2  was  380,521,  while  there  were  380,088  CSFs  of  Bj,  B2,  and  B3 
symmetry.  Compared  with  the  small-core  uranyl  MR-CISD  calculation  with  nearly  9 
million  CSFs,  the  uranate  Cl  calculation  is  quite  small.  While  this  is  not  enough  CSFs  to 
give  a  quantitative  picture  of  the  excited  states  of  uranate,  it  should  suffice  to  give  a 
qualitative  picture  of  the  nature  and  a  rough  ordering  of  excited  states.  Table  18  lists  the 
double  group  state  assignments  and  contributing  TS  terms  for  the  states  of  A  symmetry  in 
D2,  while  Table  19  lists  the  double  group  state  assignments  and  contributing  terms  for  the 
states  of  Bi  symmetry  in  D2.  Table  20  lists  the  lowest  seven  double-group  states,  their 
contributing  TS  terms,  and  the  vertical  electronic  transition  energies  to  the  ground-state. 
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Table  18.  Small-core  Uranate  MR-CISD  Double  Group  Terms  Of  A  Symmetry  in  D2 


and  Their  TS  Compositions 


Double  Group  Term 

Small-Core  Wave  Function  TS 
Composition 

1A, 

96%  'A,  from  (5t2)6 

IE 

85%  3T2  from  (5t2)5(4a1)1 

1A2 

91%  3T2  from  (5t2)5(6t2)‘ 

2A2 

55%  %  from  (5t2)5(6t2)' 

33%  3Ti  from  (5t2)5(6t2)‘ 

2E 

66%  3Ti  from  (5t2)5(6t2)‘ 

21%  JE  from  (5t2)5(6t2)1 

Table  19.  Small-core  Uranate  MR-CISD  Double  Group  Terms  Of  Bi  Symmetry  in  D2 

and  Their  TS  Compositions 


Double  Group  Term 

Small-Core  Wave  Function  TS 
Composition 

Hi 

37%  3E  from  (5t2)'(6t2)1 

26%  3T2  from  (SU)5^)1 

13%  3Ti  from  (5t2)5(6t2)‘ 

13%  3T2  from  (5t2)5(6t2)‘ 

H2 

30%  3E  from  (5t2)5(6t2)1 

28%  3Ti  from  (5t2)5(6t2)‘ 

21%  3T2  from  (5t2)5(6t2)1 

2Ti 

45%  3E  from  (5t2)5(6t2)‘ 

36%  1T1  from  (5t2)5(6t2)‘ 

3Tj 

56%  3T !  from  (5t2)5(4ai)1 

22%  3Ai  from  (5t2)5(6t2)1 

2T2 

39%  3T2  from  (5t2)5(6t2)‘ 

29%  3T !  from  (5t2)5(6t2)‘ 

19%  'T2  from  (5t2)5(6t2)‘ 

3T2 

40%  'T2  from  (SU)5^)1 

31%  3T2  from  (5t2)5(4ai)1 
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Table  20.  Small-core  Uranate  MR-CISD  Low-Lying  Vertical  Electronic  Transition 

Energies  at  1.9253  A 


Double  Group  Term 

Contributing 
TS  States 

Vertical  AE  (cm' 

2E 

66%  3T[ 
21%  lE 

34043 

2Ai 

96%  'A, 

33959 

ia2 

91%  3T2 

33815 

IE 

85%  3T2 

33571 

2Ti 

45%  3E 
36%  'T| 

33480 

it2 

30%  3E 

28%  3Ti 

21%  3T2 

33305 

1 T  i 

37%  3E 
26%  3T2 
13%  3Ti 

13%  3T2 

32633 

1A, 

96%  'A! 

0 
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V.  Conclusions  and  Recommendations 


This  thesis  has  outlined  the  theory  behind  several  theoretical  methods  used  in 
calculating  molecular  properties  of  actinide  compounds.  Excited  state  spectra  of  the 
uranyl  ion  was  calculated  using  both  large-  and  small-core  RECPs,  and  the  results  of 
these  calculations  were  compared  with  previous  theoretical  work  and  experimental 
measurements  of  fluorescent  series  spectra. 

The  large-core  uranyl  calculations  perfonned  in  this  research  faithfully 
reproduced  the  calculations  performed  at  OSU  (Zhang,  1999),  with  some  small 
variations.  These  variations  are  explained  by  the  slight  difference  in  the  treatment  of  the 
ground  and  excited  states  in  this  calculation  as  compared  with  the  OSU  calculation. 

When  compared  with  the  precise  spectroscopy  obtained  for  CS2UO2CI4  by 
Denning  et  al  (Denning  et  al,  1976;  Denning,  1992),  the  small-core  RECP  and  basis  set 
yield  a  better  match  with  experiment  on  average  than  the  large-core  results,  although  both 
RECPs  and  basis  sets  compare  favorably  with  experiment.  Because  of  this,  the  small- 
core  RECP  outperforms  the  large-core  RECP  when  examining  fluorescent  spectra  of  the 
uranyl  ion. 

The  match  between  experimental  and  theoretical  electronic  spectra  for  both  the 
large-  and  small-core  calculations  for  uranyl  validates  the  assumed  valence-core 
separability  inherent  in  the  relativistic  effective  core  potential  approximation.  The  large- 
core  uranyl  calculation  compared  reasonably  well  with  experimental  measurements  of 
excited  electronic  spectra,  and  the  calculations  involved  a  relatively  modest 
computational  effort.  The  small-core  calculation  removed  the  uranium  5d  shell  electrons 
from  the  core  potential,  and  allowed  them  to  participate  in  the  MR-CISD  calculation. 
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Assuming  the  crystal  field  perturbations  to  the  uranyl  electronic  energies  over  the  range 
of  uranyl-containing  crystals  measured  average  to  zero,  the  small  core  theoretical  result  is 
in  good  agreement  with  the  fluorescent  series  measurements  averaged  over  all  the 
experiments.  The  implication  is  then  that  the  average  experimental  fluorescent  spectra 
represents  an  unperturbed  uranyl  ion,  at  least  to  first  order.  Thus,  the  main  trade-off 
between  the  large-  and  small-core  calculation  is  the  availability  of  more  electrons  to  be 
correlated,  and  that  this  extra  electron  correlation  gives  improved  results  when  compared 
with  experiment.  Thus,  in  the  uranyl  ion,  electron  correlation  is  more  important  than 
relativistic  effects,  and  the  5d  electron  shell  in  uranium  plays  a  small  role  in  the  chemical 
bonding  in  uranium  oxides. 

The  next  reasonable  step  is  to  examine  the  excited  state  spectra  of  the  uranate 
(UO4  ')  ion  using  this  small-core  RECP.  This  ion  is  also  a  closed-shell  molecule,  and 
geometry  optimizations  using  DFT  and  MP2  indicate  the  ion  can  be  stable  in  a 
tetrahedral  geometry,  but  the  flattened  tetrahedral  geometry  had  a  slightly  lower  energy 
in  the  DFT  optimization.  Using  tetrahedral  geometry,  the  potential  energy  surfaces  could 
be  parameterized  by  a  single  parameter,  the  U-0  bond  length,  greatly  simplifying  the 
study  of  the  potential  energy  surfaces.  This  single  bond-length  represents  the  symmetric 
stretch  nonnal  mode  coordinate  of  the  molecule.  Varying  this  coordinate  simultaneously 
varies  all  four  oxygen-uranium  bond  lengths  by  an  equal  amount.  Comparison  between 
the  low-lying  excited  states  of  uranyl  and  uranate  illustrates  the  effect  the  oxygen 
coordination  has  on  the  electronic  spectra  of  uranium  oxides. 

In  the  uranyl  ion,  there  are  non-bonding  uranium  5f  orbitals,  which  are  relatively 
unaffected  by  the  oxygen  ligands.  In  the  tetrahedral  uranate  ion,  this  is  not  the  case. 

Here,  there  orbitals  are  either  bonding  or  anti-bonding:  no  non-bonding  uranium  5f 
orbitals  are  present.  The  preliminary  result  of  this  is  that  low-lying  uranate  electronic 
transitions  occur  in  the  ultraviolet  spectrum,  while  the  similar  transitions  in  the  uranyl  ion 
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are  in  the  optical  spectrum.  This  is  the  reason  for  the  characteristic  green-yellow 
fluorescence  in  many  uranium  oxides.  The  presence  or  absence  of  the  stable  uranyl  unit, 
and  thus  the  non-bonding  uranium  5f  orbitals,  dictates  the  nature  of  the  electronic 
spectra. 

The  results  of  this  research  reveal  the  methods  and  challenges  behind  the 
theoretical  study  of  actinide  molecules.  The  study  of  simple,  heavy-element  molecules 
such  as  the  uranyl  ion  has  been  quite  successful;  however,  application  of  these  theoretical 
techniques  to  more  complicated,  and  more  relevant  uranium  oxide  compounds  is 
difficult,  both  computationally,  and  theoretically.  Advances  in  parallel  computing  and 
software  resources  has  begun  to  allow  the  accurate  study  of  increasingly  complicated 
actinide  molecules.  This  capability  will  have  enormous  impact  on  such  fields  as  nuclear 
reactor  fuel  design,  nuclear  stockpile  stewardship,  and  nuclear  forensics. 


Appendix  A:  Symmetry  Considerations  for  the  Linear  Uranyl  Ion 


All  uranyl  geometry  optimization  calculations  performed  using  NWChem  4.0.1 
resulted  in  a  linear  molecule.  The  most  appropriate  symmetry  point  group  for  a  linear 
molecule  is  the  Da*  point  group.  COLUMBUS  calculations  using  the  graphical  unitary 
group  approach  (GUGA)  Cl  program  must  be  performed  using  an  abelian  point  group. 
The  abelian  sub  group  of  Dah  used  in  this  research  was  the  D2h  point-group.  The 
correlation  table  between  the  Da*  and  D2h  symmetry  point  groups  is  given  in  Table  21. 

Next,  the  symmetries  of  the  uranium  atomic  orbitals  in  Dooh  and  D2h,  as  well  as  the 
symmetries  of  linear  combinations  of  the  two  oxygen  atomic  orbitals  in  both  the  were 
found  Da)h  and  D2h  point-groups  were  found  in  order  to  create  the  ARGOS  input  fde. 
Table  22  and  Table  23  show  these  orbital  symmetries. 

Table  2 1 .  Correlation  Between  Dooh  and  D2h  Symmetry  Point-Groups  (Cotton,  1971 : 


359-362). 


Dcoh 

D2I1 

y  + 

> 

OQ 

Ag 

Ag+Big 

0g 

B2g+B3g 

89 


Table  23.  Two  Oxygen  Atomic  Orbitals  Along  z-axis  in  Dooh  and  D2h 


Orbital 

Dcoh 

P2h 

5  +  5 

+ 

°g 

ag 

5-5 

+ 

au 

biu 

Pz  +P= 

+ 

au 

biu 

Pz  ~Pz 

+ 

ag 

Py  +  Py 

Px  +  Pz 

7lu 

b2u+b3u 

Py  -Py 

Px  -  Px 

Tig 

b2g+b3g 

Because  the  references  used  in  the  MR-CISD  calculation  are  all  products  of 
ungerade  orbitals,  all  the  possible  states  arising  from  these  references  will  be  gerade. 
The  references  used  in  the  MR-CISD  calculations  were  the  3oy  ,  3 cr  ’  1  Slu ,  and  3cr 1  (p\ 

electron  configurations.  Table  24  shows  the  possible  states  arising  from  these  electron 
configurations  in  both  AS9  and  coco10  coupling  schemes. 


s  ungerade  =  odd  parity 

9  Russell-Sanders  coupling  is  called  LS  coupling  in  atoms,  AS  coupling  in  linear  molecules 

10  j-j  coupling  in  atoms,  coco  coupling  in  linear  atoms 
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Table  24.  Uranyl  Possible  States  from  CISD  References  in  both  AS  and  coco  Coupling 


Schemes 


Reference 

confieuratio 

n 

Possible  States 
(AS  coupline) 

Possible  States 
(  coco  counline) 

3<t„2 

A; 

3^:  is'H 

^gX 

Xg 

3g 

^4g 

2g 

3  <rlM 

3®3g>3, 

3g 

^g 

4g 
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Table  25  shows  the  spatial,  spin,  and  total  wave  function  symmetries  arising  from 


each  uranyl  reference  configuration. 


Table  25.  Total  Wave  Function  Symmetry  in  D2h  From  Reference  Electronic 

Configurations 


Electron 

AS 

D7h  Spatial 

D7h  Spin- 

D7h  Double-Group 

configuration 

Term 

Symmetry 

■SHmirraraB 

Symmetry 

3gu2 

+ 

Ag 

As 

Ag 

Ag 

3a*1 18,,1 

3Ag 

Ag+Big 

Big+B2g+B3g 

Big 

2  B2g 

2  B3g 

2  Ag 

So^lcpu1 

3®g 

B2g+B3g 

Big+B2g+B3g 

2  Big 

B2g 

_ B3g _ 
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Appendix  B:  Symmetry  Considerations  of  the  Tetrahedral  Uranate  Ion 


The  optimum  uranate  ion  geometry  found  using  NWChem  4.0.1  was  dependent 
on  the  starting  geometry.  If  the  geometry  optimization  was  started  in  a  tetrahedral 
configuration,  the  final  geometry  was  tetrahedral.  The  geometry  obtained  from  an 
asymmetric  starting  point  was  a  flattened  tetrahedral,  or  a  D2d  molecule.  The  D2d 
geometry  had  was  at  a  slightly  lower  energy  than  the  tetrahedral  geometry  using  DFT. 
Analysis  of  the  potential  energy  surfaces  of  a  D2d  molecule  is  complicated  by  the  fact  that 
there  are  two  independent  bond  lengths.  For  this  reason,  the  tetrahedral  molecular 
geometry  was  examined.  Using  a  tetrahedral  symmetry,  the  potential  energy  surfaces 
could  be  characterized  by  a  single  parameter,  as  was  the  case  in  the  uranyl  ion. 

A  tetrahedral  geometry  is  best  described  using  a  Td  symmetry  point  group. 
However,  as  was  the  case  with  the  uranyl  calculation,  the  COLUMBUS  calculations 
using  the  GUGA  Cl  program  must  be  perfonned  using  an  abelian  point  group.  The 
uranate  calculations  were  perfonned  using  the  D2  symmetry  point-group.  The  correlation 
table  between  the  Td  and  D2  symmetry  point  groups  is  given  in  Table  26. 
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Table  26.  Correlation  Table  Between  Td  and  D2  Symmetry  Point-Groups  (Cotton,  1971 : 


356-363). 


Id 

\h 

Ai 

A 

a2 

A 

E 

A+A 

Ti 

Bi+B2+B3 

t2 

Bi+B2+B3 

In  a  similar  fashion  to  the  uranyl  calculations,  the  Td  and  D2  symmetries  of  the 
uranium  and  four  oxygen  atomic  orbitals  was  found,  and  the  appropriate  ARGOS  input 
fde  was  generated.  The  symmetries  of  the  uranium  orbitals  in  Td  and  D2  symmetry  point- 
groups  are  shown  in  Table  27,  while  the  symmetries  of  the  linear  combinations  of  four 
oxygen  atoms  are  shown  in  Table  28. 
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Table  28.  Combination  of  Four  Tetrahedral  Oxygen  Atomic  s-Orbitals  in  T<i  and  D2 


Oxveen  orbital 

Id 

D2 

s 

Ai+T2+T! 

A+B 1+B2+B3 

P: 

Ai+T2+T! 

A+B 1+B2+B3 

Py 

Ai+T2+T! 

A+B  1+B2+B3 

Px 

A!+T2+T! 

A+B1+B2+B3 

Three  electronic  configurations  were  used  as  references  in  the  MR-CISD 
calculation:  the  5 1\  ground  state,  as  well  as  5t\%a\  and  5 t\6t\ ,  which  represent 
excitations  in  the  uranium  5f  orbitals.  The  possible  double  group  states  in  Td  and  D2 
symmetry  point  groups  arising  from  these  three  references  are  listed  in  Table  29 
(Herzberg,  1966:  570). 
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Table  29.  Uranate  Total  Wave  Function  Symmetry  in  D2  From  Reference  Electronic 


Configurations 


Electron 

configuration 

Td  rs11 

Terms 

Tri  Spatial 
Symmetry 

Tri  Spin- 
Symmetry 

Td  Double 
Group  Terms 

D?  Double  Group 
Terms 

5*26 

lA 

Ai 

Ai 

Ai 

A 

5^  8a} 

% 

t2 

Ai 

t2 

B1+BAB3 

a2 

A 

E 

2A 

5*2  8a' 

% 

t2 

Ti 

Tt 

B1+BAB3 

t2 

Bi+B2+B3 

% 

Ai 

A, 

A 

5*2  6*2 

'E 

E 

E 

A+A 

Ai 

't 

1\ 

Ti 

Tt 

B1+BAB3 

% 

t2 

t2 

B1+BAB3 

Ai 

Ti 

B1+BAB3 

3e 

E 

Ti+T2 

2Bi+2B2+2B3 

5*2  6*2 

o _ 

Ti 

3t 

1 1 

Ti 

A!+E+Ti+T2 

3A+2B  i+2B2+2B3 

3t 

12 

t2 

A2+E+Ti+T2 

3A+2B  i+2B2+2B3 

Total 

3  A  i+2A2+4E+ 

6T1+6T1 

1 3  A+ 1 2(B  i+B2+B3) 

It  is  evident  from  Table  29  that  the  total  number  of  possible  states  arising  from  the 
three  electron  configuration  references  is  much  larger  in  the  tetrahedral  uranate  ion  than 
in  the  linear  uranyl  ion. 


n 


Russell-Sanders  coupling  in  non-linear  poly  atomic  atoms  is  called  FS  coupling 
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Appendix  C:  State  Assignments  Using  MR-CISD  Wave  Functions 


For  both  the  uranyl  and  uranate  ions,  the  MR-CISD  calculation  was  performed 
using  an  abelian  point-group  that  was  a  sub-group  of  the  actual  symmetry  group 
corresponding  to  the  molecular  symmetry.  The  SCF,  MCSCF  and  MR-CISD  wave 
functions  obtained  in  the  calculations  all  reflect  the  symmetry  of  the  calculation,  which  in 
turn,  echoed  the  overall  molecular  symmetry. 

In  the  MR-CISD  calculations,  symmetry  considerations  and  degeneracy  reduced 
the  total  number  of  calculations  necessary  to  characterized  all  the  possible  states  arising 
from  the  electron  references.  For  the  uranyl  calculation,  calculating  four  MR-CISD  roots 
of  Ag  symmetry  and  three  roots  of  B2g  symmetry  completely  characterized  the  ground 
state  and  six  excited  electronic  states  possible  from  reference  electron  configurations  (see 
Table  24  and  Table  25).  Once  these  MR-CISD  wave  functions  are  obtained,  all  that 
remains  is  identification  of  each  state  corresponding  to  the  particular  MR-CISD  root. 

Uranyl  MR-CISD  State  Assignment 

In  the  case  of  the  uranyl  ion,  the  MR-CISD  state  assignment  was  relatively 
straightforward,  using  the  infonnation  in  Table  24  and  the  process  of  elimination. 

The  output  of  the  MR-CISD  calculation  lists  electron  configurations  in  the  Cl 
expansion  with  a  Cl  coefficient  above  a  certain  threshold.  Also  listed  is  the  spin 
multiplicity  of  the  configuration.  Listed  below  is  an  excerpt  from  the  first  MR-CISD  root 
of  Ag  symmetry  in  the  small-core  uranyl  calculation. 
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—  list  of  ci  coefficients  (  ctol  =  1.00E-02  )  total  energy(  1)  =  -246.4005545447 


internal  orbitals 


level  1  2  3  4  5  6  7  8  9  10  11  12  13  14  15  16 

orbital  3  4  5  47  61  77  78  93  94  106  107  79  74  80  95  108 
symmetry  ag  ag  ag  b2g  b3g  blu  blu  b2u  b2u  b3u  b3u  blu  au  blu  b2u  b3u 


path  s  ms  csf '#  c(i)  ext.  orb.(sym) 


z* 

1 

1 

1  0.908837 

+- 

+- 

+- 

+- 

+-  +- 

+- 

+-  +- 

+-  +- 

z 

1 

1 

6  -0.035430 

+- 

+- 

+- 

+- 

+-  +- 

+- 

+-  +- 

+-  +- 

z 

1 

1 

10  -0.035430 

+- 

+- 

+- 

+- 

+-  +- 

+- 

+-  +- 

+-  +- 

z 

1 

1 

13  -0.011775 

+- 

+- 

+- 

+- 

+-  +- 

+- 

+-  +- 

+-  +- 

The  orbital  number  lists  the  unfrozen,  occupied  core  orbitals  and  reference 
orbitals.  The  next  line  lists  the  symmetry  of  each  orbital.  The  pluses  and  minuses 
represent  electron  spin-up  or  spin-down.  On  the  left  side,  the  c(i)  represents  the  CI 
coefficient,  and  the  s  represents  the  spin-multiplicity.  For  this  example,  the  79  b3u  is  the 
3cr2  HOMO  and  is  a  singlet.  This  corresponds  to  the  HF  SCF  ground  state,  and  is  the 

0 „  state.  Identification  of  the  other  states  proceeds  in  a  similar  fashion.  States  are  most 
easily  identified  through  their  Q  values  and  spin-multiplicity.  For  example,  roots  of  Ag 
symmetry  can  only  correspond  to  even  Q  values.  Thus,  Ag  roots  can  only  be  0^ ,  2g ,  or 


4g  states,  while  Big  roots  can  only  be  1  or  3„  states.  The  4a  state  can  be  identified  by 
the  fact  that  there  can  be  no  contributions  from  the  3cr'  1 8lu  configuration.  It  can  only 
arise  from  the  3cr'  \m]  reference  or  from  excitations  from  n ,  orbitals.  The  2„  and  3„ 

U  f  u  u  g  g 

states  can  be  differentiated  by  the  singlet  states.  For  example,  the  2g  states  contains  a 
small  contribution  from  the  singlet  3cr'(  1 8lu  configuration.  This  uniquely  identifies  the 
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2  states.  The  3  states  are  identified  by  the  presence  of  the  singlet  3  c,'  1  (p\ 
configuration.  The  1  states  are  identified  by  the  absence  of  any  singlet  configurations 
from  the  3  c,'  1 8\  or  3 c,1,  1  cp\  reference  configurations.  Thus,  all  seven  states  arising  from 
the  three  references  are  uniquely  identified  and  ordered  in  energy. 

Uranate  MR-CISD  State  Assignment 

State  identification  in  the  tetrahedral  uranate  ion  proceeds  in  a  similar  fashion. 

The  E  states  are  identified  by  their  degeneracy.  Two  consecutive,  degenerate  roots  of  A 
symmetry  correspond  to  an  E  state.  The  Ai  state  can  be  identified  by  noting  that  the  first 
root  of  A  symmetry  corresponds  to  the  ground  state,  a  1Al  state.  Identification  of  A2 
states  is  slightly  more  complicated,  and  it  relies  on  the  Cl  coefficients  of  the  ground  state. 
The  Ai  and  Ao  states,  while  both  of  A  symmetry  in  D2,  are  differentiable  by  the  relative 
signs  of  the  Cl  coefficients  of  various  electron  configurations  in  each  wave  function.  The 
key  to  identification  of  an  A2  state  is  to  find  electron  configurations  in  the  ground  state 
that  are  identical  in  the  A2  state,  and  examine  the  Cl  coefficients.  The  Cl  coefficients  of 
the  electron  configurations  in  the  Ai  state  should  come  into  the  A2  state  with  different 
signs,  while  they  will  have  the  same  sign  in  another  Ai  state.  Thus,  Ai  and  A2  states  can 
be  differentiated  by  the  phases  of  their  wave  functions,  relative  to  the  Cl  coefficients  in 
the  ground  state,  once  a  series  of  identical  electron  configurations  are  located  between  the 
ground  state  (Ai)  and  the  state  to  be  identified. 

Differentiation  between  Ti  and  T2  states  proceeds  along  a  similar  fashion. 
However,  the  complication  here  is  to  first  identify  either  a  Ti  or  T2  state,  then  use  the 
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relative  phases  of  the  Cl  coefficients  in  identical  electron  configurations  to  distinguish 
between  the  two.  Here,  the  spin  multiplicity  can  help.  Ti  or  T?  states  can  be  identified 
by  locating  a  singlet  state  from  a  t\a\  or  t\a\  electron  configuration.  Therefore,  Ti  and 
T2  states  can  be  identified  from  either  singlet  states  arising  from  excitations  into  ai  from 
ti  or  t2  orbitals,  or  into  fi  or  t2  orbitals  from  orbitals.  Once  one  such  state  is  found, 
other  Ti  or  T2  states  can  be  differentiated  by  examining  the  relative  phases  of  the  Cl 
coefficients. 
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List  of  Abbreviations 


AFIT 

Air  Force  Institute  of  Technology 

AMD 

Advanced  Micro  Devices 

ARECP 

Averaged  Relativistic  Effective  Core  Potential 

ASC 

Aeronautical  System  Center 

aug-cc-pvdz 

Augmented,  Correlation-Consistent  Polarized  Valence  Double  Zeta 

B3LYP 

Becke  Three  Parameter  Lee- Yang-Parr  hybrid  exchange-correlation 

functional 

CAS 

Complete  Active  Space 

CCSD(T) 

Coupled-Cluster  Singles  and  Doubles  (Triples) 

CISD 

Configuration  Interaction,  Single  and  Double  excitations 

CSF 

Configuration  State  Function 

DFT 

Density  Functional  Theory 

DHF 

Dirac-Hartree-Fock 

EMSL 

Environmental  Molecular  Sciences  Laboratory 

eV 

electron  volt  ( 1 . 602 1 76  x  1 0' 1 9  J) 

Gb 

Gigabyte 

GCA 

Gradient-Corrected  Approximation 

GHz 

Gigahertz 

HF 

Hartree-Fock 

HOMO 

Highest  Occupied  Molecular  Orbital 

kcal 

kilocalories  (4184  J) 

LUMO 

Lowest  Unoccupied  Molecular  Orbital 

Mb 

Megabyte 
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MCSCF 


MHz 

MP2 

MR-CISD 

MSRC 

O 

OSU 

PES 

PC 

PNNL 

RAM 

RECP 

SCF 

SO 

SOCI 

Tb 

U 


Multi-Configuration  Self-Consistent  Field 
Megahertz 

Moller-Plesset  second-order  perturbation  theory 
Multi-Reference  Configuration  Interaction  with  Single  and  Double 
excitations 

Major  Shared  Resource  Center 
Oxygen 

Ohio  State  University 
Potential  Energy  Surface 
Personal  Computer 

Pacific  Northwest  National  Laboratory 
Random  Access  Memory 
Relativistic  Effective  Core  Potential 
Self-Consistent  Field  method 
Spin-Orbit 

Spin-Orbit  Configuration  Interaction 

Terabyte 

Uranium 
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